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Abstract. We report some recent developments in the understanding of GRBs based on the theoretical framework of the 
"fireshell" model, already presented in the last three editions of the "Brazilian School of Cosmology and Gravitation". After 
recalling the basic features of the "fireshell model", we emphasize the following novel results: 1) the interpretation of the 
X-ray flares in GRB afterglows as due to the interaction of the optically thin fireshell with isolated clouds in the CircumBurst 
Medium (CBM); 2) an interpretation as "fake - disguised" short GRBs of the GRBs belonging to the class identified by Norris 
& Bonnell; we present two prototypes, GRB 970228 and GRB 060614; both these cases are consistent with an origin from 
the final coalescence of a binary system in the halo of their host galaxies with particularly low CBM density n c i, m ~ 
particles/cm 3 ; 3) the first attempt to study a genuine short GRB with the analysis of GRB 050509B, that reveals indeed still 
an open question; 4) the interpretation of the GRB-SN association in the case of GRB 060218 via the "induced gravitational 
collapse" process; 5) a first attempt to understand the nature of the "Amati relation", a phenomenological correlation between 
the isotropic-equivalent radiated energy of the prompt emission E[ so with the cosmological rest-frame vFy spectrum peak 
energy E p j. In addition, recent progress on the thermalization of the electron-positron plasma close to their formation phase, as 
well as the structure of the electrodynamics of Kerr-Newman Black Holes are presented. An outlook for possible explanation 
of high-energy phenomena in GRBs to be expected from the AGILE and the Fermi satellites are discussed. As an example 
of high energy process, the work by Enrico Fermi dealing with ultrarelativistic collisions is examined. It is clear that all the 
GRB physics points to the existence of overcritical electrodynamical fields. In this sense we present some progresses on a 
unified approach to heavy nuclei and neutron stars cores, which leads to the existence of overcritical fields under the neutron 
star crust. 



INTRODUCTION 

Gamma-Ray Bursts (GRBs) represent very likely "the" most extensive computational, theoretical and observational 
effort ever carried out successfully in physics and astrophysics. The extensive campaign of observation from space 
based X-ray and y-ray observatory, such as the Vela, CGRO, BeppoS AX, HETE-II, INTEGRAL, Swift, Agile, GLAST, 
R-XTE, Chandra, XMM satellites, have been matched by complementary observations in the radio wavelength (e.g. 
by the VLA) and in the optical band (e.g. by VLT, Keck, REM). The very fortunate situation occurs that these data can 
be confronted with a mature theoretical development. 

We outline how this progress leads to the confirmation of three interpretation paradigms for GRBs we proposed 



1 Part I, Part II and Part III of these Lecture notes have been published respectively in COSMOLOGY AND GRAVITATION: X"' Brazilian School of 
Cosmology and Gravitation; 25 th Anniversary (1977-2002), M. Novello, S.E. Perez Bergliaffa (eds.), AIP Conf. Proc, 668, 16 (2003), see Ref. [1], 
in COSMOLOGY AND GRAVITATION: XT' 1 Brazilian School of Cosmology and Gravitation, M. Novello, S.E. Perez Bergliaffa (eds.), AIP Conf. 
Proc, 782, 42 (2005), see Ref. [2], and in COSMOLOGY AND GRAVITATION: XII th Brazilian School of Cosmology and Gravitation, M. Novello, 
S.E. Perez Bergliaffa (eds.), AIP Conf. Proc, 910, 55 (2007), see Ref. [3]. 



in July 2001 [4, 5, 6]. The outcome of this analysis points to the existence of a "canonical" GRB, originating from 
a variety of different initial astrophysical scenarios. The communality of these GRBs appears to be that they all are 
emitted in the process of formation of a black hole with a negligible value of its angular momentum. The following 
sequence of events appears to be canonical: the gravitational collapse to a black hole, the vacuum polarization process 
in the dyadosphere with the creation of the optically thick self accelerating electron-positron plasma; the engulfment of 
baryonic mass during the plasma expansion; adiabatic expansion of the optically thick "fireshell" of electron-positron- 
baryon plasma up to the transparency; the interaction of the accelerated baryonic matter with the CircumBurst Medium 
(CBM). This leads to the canonical GRB composed of a proper GRB (P-GRB), emitted at the moment of transparency, 
followed by an extended afterglow. The sole parameters in this scenario are the total energy of the dyadosphere E^y a , 
the fireshell baryon loading M B defined by the dimensionless parameter B = M B c 2 /E dya , and the CBM filamentary 
distribution around the source. In the limit B —> the total energy is radiated in the P-GRB with a vanishing contribution 
in the extended afterglow. We refer globally to this model as the "fireshell" model. 

The increase of observational details on many different GRBs, gained in the last two year observational campaign, is 
showing very clearly an evolution in the understanding of this basic scenario. This has allowed to explore GRBs with 
different baryon loading and very different CBM properties. A major new result has been obtained in the understanding 
of the origin of flares, which is traced back to the interaction of the fireshell with isolated CBM clouds (see below, 
section GRB 060607 'A: a complete analysis of the prompt emission and X-ray flares). 

Another key result, presented in these lectures, has been the understanding of a new GRB class, which was 
pioneered by Norris and Bonnell [7], characterized by "an occasional softer extended emission lasting tenths of 
seconds after an initial spikelike emission". From the "fireshell" model it clearly follows that these sources explore 
a new range of parameters and mark a fundamental step in the understanding of the astrophysical nature of the GRB 
progenitor systems. The new class is shown to be characterized by GRBs occurring in a particularly low density CBM, 
n cbm ~ 10~ 3 particles/cm 3 , typical of galactic halos. The progenitors must therefore necessarily be binary sources 
which have migrated from their birth location in a star forming region to a low density region within the galactic halo, 
where the final merging occurs. If such sources did explode in a CBM with average density n c ), m <~ 1 particle/cm 3 , they 
would look exactly like what is called long GRBs. In this sense, we decided to call them "fake" or "disguised" short 
GRBs (see below, section The Norris & Bonnel kind of sources: the new class of "fake - disguised" short GRBs). 

This leads, in turn, to a novel interpretation of the traditional GRB classification [8, 9] in "short" GRBs (with a Tgo 
duration lasting less than ~2 s) and "long" ones (with a 790 duration lasting more than <~ 2 s up to <~ 1000 s). The 
duration of long GRBs is shown to be just due to the instrumental noise threshold and not to represent any intrinsic 
GRB feature. On the other hand, there is an increasing evidence that the majority of the sources classified as short 
GRBs are actually "fake - disguised" short ones. No Supernova can possibly be related to these GRBs (see below, 
section The "fireshell" model and GRB progenitors). However, some issues remains open (see below, section Open 
issues in current theoretical models). Moreover, the quest for the first clear identification of a "genuine" short GRB is 
also still open (see below, section The search for a "genuine short" GRB: the case of GRB 050509B). 

The GRBs associated with SN do necessarily form a different class, of weakest and more numerous sources, 
originating also from binary systems, formed by a neutron star, close to its critical mass, and a companion star, 
evolved out of the main sequence, via the "induced gravitational collapse" process (see below, section GRBs and 
SNe: the induced gravitational collapse). 

The detection of GRBs up to very high redshifts [up to z = 6.7, see Ref. 10], their high observed rate of one every 
few days, and the progress in the theoretical understanding of these sources all make them useful as cosmological 
tools, complementary to supernovae la, which are observed only up to z = 1.7 [11, 12]. One of the hottest topics 
on GRBs is the possible existence of empirical relations between GRB observables [13, 14, 15, 16, 17, 18], which 
may lead, if confirmed, to using GRBs as tracers of models of universe. The first empirical relation, discovered when 
analyzing the BeppoSAX so-called "long" bursts with known redshift, was the "Amati relation" [13]. It was found that 
the isotropic-equivalent radiated energy of the prompt emission Ej so is correlated with the cosmological rest-frame 
vF v spectrum peak energy E p /. E p i °< (Ei so ) a , with a = 0.52 ±0.06 [13]. The existence of the Amati relation has 
been confirmed by studying a sample of GRBs discovered by Swift, with a = 0.491qo5 [19, 20]. We present a first 
attempt to understand the nature of such a phenomenological correlation (see below, section Theoretical background 
for GRBs' empirical correlations). 

All these theoretical and phenomenological approach has motivated significant progresses on: 

a) The electron-positron plasma thermalization. We solved numerically relativistic Boltzmann equations with 
collisional integrals representing two-body as well as three-body interactions, in particular Compton, Moller and 
Bhabha scattering, pair creation and annihilation, relativistic bremsstrahlung, three photon creation and annihilation, 




double Compton scattering and the corresponding processes where protons participate. This allowed determination 
of characteristic timescales of thermalization as well as clarified the role of binary and triple interactions in reaching 
thermal equilibrium (see below, section Thermalization process of electron-positron plasma with baryon loading). 

b) The electrodynamics of neutron stars and black holes. We present a unified treatment of nuclear density cores 
recovering the classical results of neutral atoms with heavy nuclei with mass number A w 10 2 -10 6 and extrapolating 
these results to massive nuclear density cores with A rts {mpianck^nf 1 ~ 10 57 . The treatment is approached by solving 
the relativistic Thomas-Fermi equation describing a system of N n neutrons, N p protons and N e electrons in beta 
equilibrium. In order to show the stability of such cores under the competing effects of self gravity and Coulomb 
repulsion, we have started to take the gravitational field into duly account, and put the issue within the framework of 
general relativity. We expect that starting from such configuration, gravitational collapse would lead to Dyadosphere 
(electron-positron-photon plasma) in the Reissner-Nordstrm geometry, or Dyadotorus in the Kerr-Newmann geometry. 
We present in some details the analysis of the Dyadotorus (see below, section Critical electric Fields on the surface of 
massive cores and Dyadotorus of the Kerr-Newman Geometry). 

c) In view of the new data from the Fermi and AGILE satellites, an analysis of the GRB radiation over 1 MeV. It is 
by now clear that our pure thermal emission previously considered, and which has been fundamental in expressing the 
average CBM density, is not appropriate to the description of this high-energy component. In parallel, we are currently 
examining how Fermi ideas [21] have been further developed in large data analysis procedures at CERN and other 
accelerators all over the world (see below, section Selected processes originating high-energy emission). 



BRIEF REMINDER OF THE FIRESHELL MODEL 

The black hole uniqueness theorem [see Left panel in Fig. 1 and e.g. Ref. 22] is at the very ground of the fact 
that it is possible to explain the different Gamma-Ray Burst (GRB) features with a single theoretical model, over 
a range of energies spanning over 6 orders of magnitude. The fundamental point is that, independently of the fact that 
the progenitor of the gravitational collapse is represented by merging binaries composed by neutron stars and white 
dwarfs in all possible combinations, or by a single process of gravitational collapse, or by the process of "induced" 
gravitational collapse, the formed black hole is totally independent from the initial conditions and reaches the standard 
configuration of a Kerr-Newman black hole (see Right panel in Fig. 1). It is well known that pair creation by vacuum 
polarization process can occur in a Kerr-Newman black hole [23, 24]. 

We consequently assume, within the fireshell model, that all GRBs originate from an optically thick e plasma 
with total energy Ef ot in the range 10 49 -10 54 ergs and a temperature T in the range 1-4 MeV [26]. Such an e ± 
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Figure 2. The "canonical GRB" light curve theoretically computed for GRB 991216. The prompt emission observed by BATSE 
is identified with the peak of the extended afterglow, while the small precursor is identified with the P-GRB. For this source we 
have E'fJ = 4.83 x 1(P 3 ergs, B ~ 3.0 x 1(T 3 and (n cbm ) ~ 1.0 particles/cm 3 . Details in Ruffini et al. [5, 25, 3]. 



plasma has been widely adopted in the current literature [see e.g. Refs. 27, 28, and references therein]. After an early 
expansion, the e -photon plasma reaches thermal equilibrium with the engulfed baryonic matter Mb described by the 
dimensionless parameter B = Mbc 2 /Ef ot , that must be B < 10~ 2 [29, 30]. As the optically thick fireshell composed 
by e ± -photon-baryon plasma self-accelerate to ultrarelativistic velocities, it finally reaches the transparency condition. 
A flash of radiation is then emitted. This is the P-GRB [5]. Different current theoretical treatments of these early 
expansion phases of GRBs are compared and contrasted in Bianco et al. [31] and Ruffini et al. [32]. The amount of 
energy radiated in the P-GRB is only a fraction of the initial energy Ef ot . The remaining energy is stored in the kinetic 
energy of the optically thin baryonic and leptonic matter fireshell that, by inelastic collisions with the CBM, gives rise 
to a multi-wavelength emission. This is the extended afterglow. It presents three different regimes: a rising part, a peak 
and a decaying tail. We therefore define a "canonical GRB" light curve with two sharply different components (see 
Fig. 2) [5, 3, 33, 34, 35]: 1) the P-GRB and 2) the extended afterglow. What is usually called "Prompt emission" in the 
current literature mixes the P-GRB with the raising part and the peak of the extended afterglow. Such an unjustified 
mixing of these components, originating from different physical processes, leads to difficulties in the current models 
of GRBs, and can as well be responsible for some of the intrinsic scatter observed in the Amati relation [20, 36, , see 
also below, section Theoretical background for GRBs' empirical correlations]. 



The optically thick phase 

In Fig. 3 we present the evolution of the optically thick fireshell Lorentz y factor as a function of the external radius 
for 7 different values of the fireshell baryon loading B and two selected limiting values of the total energy E"± of the 
e plasma. We can identify three different eras: 

1 . Era I: The fireshell is made only of electrons, positrons and photons in thermodynamical equilibrium (the "pair- 
electromagnetic pulse", or PEM pulse for short). It self-accelerate and begins its expansion into vacuum, because 
the environment has been cleared by the black hole collapse. The Lorentz 7 factor increases with radius and the 
dynamics can be described by the energy conservation and the condition of adiabatic expansion [29, 31]: 

r 0v v = o (i) 

where T^ v is the energy-momentum tensor of the e + e~ plasma (assumed to be a perfect fluid), e is its internal 
energy density, V and "f are its volumes in the co-moving and laboratory frames respectively, T is the thermal 
index and the quantities with and without the o subscript are measured at two different times during the expansion. 

2. Era II: The fireshell impacts with the non-collapsed bayonic remnants and engulfs them. The Lorentz 7 factor 
drops. The dynamics of this era can be described by imposing energy and momentum conservation during the 
fully inelastic collision between the fireshell and the baryonic remnant. For the fireshell solution to be still valid, 
it must be B< 10~ 2 [30]. 

3. Era III: The fireshell is now made of electrons, positrons, baryons and photons in thermodynamical equilibrium 
(the "pair-electromagnetic -baryonic pulse", or PEMB pulse for short). It self-accelerate again and the Lorentz 
7 factor increases again with radius up to when the transparency condition is reached, going to an asymptotic 
value Yasym = 1/B. If B ~ 10~ 2 the transparency condition is reached when 7 ~ Jasym- On the other hand, when 
B < 10~ 2 , the transparency condition is reached much before 7 reaches its asymptotic value [30, 5, , see next 
section for details]. In this era the dynamical equations are the same of the first one, together with the baryon 
number conservation: 
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where n B is the baryon number density in the fireshell. In this era it starts to be crucial the contribution of the rate 
equation to describe the annihilation of the e + e~ pairs: 



where N e ± is the number of e+e pairs and N e ±(T) is the number of e + e pairs at thermal equilibrium at 
temperature T [30, 31]. 

In the "fireball" model in the current literature the baryons are usually considered present in the plasma since 
the beginning. In other words, in the fireball dynamics there is only one era corresponding to the Era III above 
[37, 38, 39, 31]. Moreover, the rate equation is usually neglected, and this affects the reaching of the transparency 
condition. A detailed comparison between the different approaches is reported in Bianco et al. [31]. 



At the transparency point, the value of the B parameter rules the ratio between the energetics of the P-GRB and 
the kinetic energy of the baryonic and leptonic matter giving rise to the extended afterglow. It rules as well the time 
separation between the corresponding peaks [5, 32]. 

We have recently shown [40] that a thermal spectrum still occurs in presence of e pairs and baryons. By solving 
the rate equation we have evaluated the evolution of the temperature during the fireshell expansion, all the way 
up to when the transparency condition is reached [29, 30]. In the upper panel of Fig. 4 we plot, as a function of 
B, the fireshell temperature T a at the transparency point, i.e. the temperature of the P-GRB radiation. The plot is 
drawn for four different values of E t0 ± in the interval [10 49 , 10 55 ] ergs, well encompassing GRBs' observed isotropic 
energies. We plot both the value in the co-moving frame T"" n and the one Doppler blue-shifted toward the observed 
jobs = (i + j5 )YoT"' m , where j3 D is the fireshell speed at the transparency point in units of c [30]. 

In the middle panel of Fig. 4 we plot, as a function of B, the fireshell Lorentz gamma factor at the transparency 
point y a . The plot is drawn for the same four different values of E'± of the upper panel. Also plotted is the asymptotic 
value y = 1 /B, which corresponds to the condition when the entire initial internal energy of the plasma E'J± has been 
converted into kinetic energy of the baryons [30]. We see that such an asymptotic value is approached for B — > 10~ 2 . 
We see also that, if E'J* increases, the maximum values of y are higher and they are reached for lower values of B. 

In the lower panel of Fig. 4 we plot, as a function of B, the total energy radiated at the transparency point in the 
P-GRB and the one converted into baryonic and leptonic kinetic energy and later emitted in the extended afterglow. 
The plot is drawn for the same four different values of E"± of the upper panel and middle panels. We see that for 
B < 10~ 5 the total energy emitted in the P-GRB is always larger than the one emitted in the extended afterglow. In the 
limit B -> it gives rise to a "genuine" short GRB (see also Fig. 5). On the other hand, for 3.0 x 10~ 4 < B < 10~ 2 
the total energy emitted in the P-GRB is always smaller than the one emitted in the extended afterglow. If it is not 
below the instrumental threshold and if n c b m ~ 1 particle/cm 3 (see Fig. 14), the P-GRB can be observed in this case as 
a small pulse preceding the main GRB event (which coincides with the peak of the extended afterglow), i.e. as a GRB 
"precursor" [5, 1, 35, 32, 33, 41]. 

Particularly relevant for the new era of the Agile and GLAST satellites is that for B < 10~ 3 the P-GRB emission 
has an observed temperature up to 10 3 keV or higher. This high-energy emission has been unobservable by the Swift 
satellite. 
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Figure 3. The Lorentz y factor of the expanding fireshell is plotted as a function of its external radius for 7 different values of the 



fireshell baryon loading B, ranging from B = 10 and B = 10 , and two selected limiting values of the total energy E ± of the 
e ± plasma: E' e g> = 1.17 x 10 49 ergs (upper panel) and E'Jl' = 1.47 x 10 53 ergs (lower panel). The asymptotic values y — > l/B are 
also plotted (dashed horizontal lines). The lines are plotted up to when the fireshell transparency is reached. 



10 4 



1CT r 



co 10* 



10' 



E 10" 



10" 



10"' 



! 10' 



Doppler blue-shifted toward the observer 



In the co-moving frame 



I 10 3 



10' 



0.8 



4J 0.6 



0.4 - 



0.2 



E^'s 1.44x10"- erg 
Eg'sl.SSxIO 51 erg 
Eg's 1.77x1 53 erg 
E^.'s; 1.22x10 55 erg 



Ef= 1. 44x1 49 erg 




Figure 4. At the fireshell transparency point, for 4 different values of E*±f, we plot as a function of B: (Above) The fireshell 
temperature in the co-moving frame r o eom (thicker lines) and the one Doppler blue-shifted along the line of sight toward the 
observer in the source cosmological rest frame T° bs (thinner lines); (Middle) The fireshell Lorentz gamma factor y together with 
the asymptotic value y D = 1 /B; (Below) The energy radiated in the P-GRB (thinner lines, rising when B decreases) and the one 
converted into baryonic kinetic energy and later emitted in the extended afterglow (thicker lines, rising when B increases), in units 
of£^ ? . 
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Figure 5. Here the energies emitted in the P-GRB (solid line) and in the extended afterglow (dashed line), in units of the total 
energy of the plasma, are plotted as functions of the B parameter for a typical value of E*°t ~ 10 53 erg (see lower panel of Fig. 4). 
When B < 10~ 5 , the P-GRB becomes predominant over the extended afterglow, giving rise to a "genuine" short GRB. In the figure 
are also marked the values of the B parameters corresponding to some GRBs we analyzed, all belonging to the class of long GRBs. 



The optically thin phase 

The dynamics of the optically thin fireshell of baryonic matter propagating in the CBM can be obtained from the 
relativistic conservation laws of energy and momentum [see e.g. Ref. 42]: 

dE mi = (y- \)dM chm c 2 
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dE m + dM cbm 



k dM chm = 4nm p n cbm r dr 

where 7, E ml and M are the pulse Lorentz gamma factor, internal energy and mass-energy respectively, n cbm is the 
CBM number density, m p is the proton mass, e is the emitted fraction of the energy developed in the collision with the 
CBM andM c b m is the amount of CBM mass swept up within the radius r. M cbm = m p n cbm (A7l /3)(r 3 — r 3 ), where r a 
is the starting radius of the shock front. 

In both our approach and in the other ones in the current literature [see e.g. Refs. 43, 44, 42, 45, 3] such conservations 
laws are used. The main difference is that in the current literature an ultra-relativistic approximation, following the 
Blandford and McKee [46] self-similar solution, is widely adopted, leading to a simple constant-index power-law 
relations between the Lorentz 7 factor of the optically thin "fireshell" and its radius: 
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with a — 3 in the fully radiative case and a — 3/2 in the adiabatic case [43, 45]. On the contrary, we use the exact 
solutions of the equations of motion of the fireshell [47, 42, 45, 48, 3]. In the adiabatic regime (e = 0) we get: 



o = jj + 27o (M cbm /M B ) + (M cbm /M B ) 2 
1 + 2y D (M cbm /M B ) + {M cbm /M B ) 2 ' 

where y D and M B are respectively the values of the Lorentz gamma factor and of the mass of the accelerated baryons 
at the beginning of the afterglow phase. In the fully radiative regime (e = 1), instead, we have: 

= 1 + (M cbm /M B ) (1 + fc*) [1 + (1/2) (M cbm /M B )} 
7 Y° l + (M cbm /M B ) (1 + Yo l ) [1 + (1/2) (M cbm /M B )} ' 

A detailed comparison between the equations used in the two approaches has been presented by Bianco and Ruffini 
[47, 42, 45, 48]. In particular, Bianco and Ruffini [45] show that the regime represented in Eq.(8) is reached only 
asymptotically when y, ^> y ^> 1 in the fully radiative regime and y 2 3> y 2 3> 1 in the adiabatic regime, where y, the 
initial Lorentz gamma factor of the optically thin fireshell. 

In Fig. 6 we show the differences between the two approaches. In the upper panel there are plotted the exact solutions 
for the fireshell dynamics in the fully radiative and adiabatic cases. In the lower panel we plot the corresponding 
"effective" power-law index a e /f, defined as the index of the power-law tangent to the exact solution [45]: 

Such an "effective" power-law index of the exact solution smoothly varies from to a maximum value which is always 
smaller than 3 or 3/2, in the fully radiative and adiabatic cases respectively, and finally decreases back to (see Fig. 6). 



Extended afterglow luminosity and spectrum 

The extended afterglow luminosity in the different energy bands is governed by two quantities associated to the 
environment. Within the fireshell model, these are the effective CBM density profile, n c \, m , and the ratio between the 
effective emitting area A e ff and the total area A tot of the expanding baryonic and leptonic shell, & = A e ff /A tot . This 
last parameter takes into account the CBM filamentary structure [49, 50] and the possible occurrence of a fragmentation 
in the shell [51]. 

Within the "fireshell" model, in addition to the determination of the baryon loading, it is therefore possible to infer 
a detailed description of the CBM, its average density and its porosity and filamentary structure, all the way from the 
black hole horizon to distance r < 10 17 cm. This corresponds to the prompt emission. This description is lacking in 
the traditional model based on the synchrotron emission. The attempt to use the internal shock model for the prompt 
emission [see e.g. Refs. 52, 27, 28, and references therein] only applies to regions where r > 10 17 cm [53]. 

In our hypothesis, the emission from the baryonic and leptonic matter shell is spherically symmetric. This allows us 
to assume, in a first approximation, a modeling of thin spherical shells for the CBM distribution and consequently to 
consider just its radial dependence [25]. For simplicity, and in order to have an estimate of the energetics, the emission 
process is postulated to be thermal in the co-moving frame of the shell [49]. We are currently examining a departure 
from this basic mechanism by taking into account inverse Compton effects due to the electron collisions with the 
thermal photons. The observed GRB non-thermal spectral shape is due to the convolution of an infinite number of 
thermal spectra with different temperatures and different Lorentz and Doppler factors. Such a convolution is to be 
performed over the surfaces of constant arrival time of the photons at the detector [EQuiTemporal Surfaces, EQTSs; 
see e.g. Ref. 42] encompassing the whole observation time [54]. 

In Fig. 7 we plot, as a function of B, the arrival time separation At a between the P-GRB and the peak of the extended 
afterglow measured in the cosmological rest frame of the source. Such a time separation At a is the "quiescent time" 
between the precursor (i.e. the P-GRB) and the main GRB event (i.e. the peak of the extended afterglow). The plot 
is drawn for the same four different values of E tc ± of Fig. 4. The arrival time of the peak of the extended afterglow 
emission depends on the detailed profile of the CBM density. In this plot it has been assumed a constant CBM density 
n c bm = 1.0 particles/cm 3 . We can see that, for 3.0 x 10~ 4 < B < 10~ 2 , which is the condition for P-GRBs to be 
"precursors" (see above), At a increases both with B and with E'J±. We can have At a > 10 2 s and, in some extreme cases 
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Figure 6. In the upper panel, the analytic behavior of the Lorentz y factor during the afterglow era is plotted versus the radial 
coordinate of the expanding optically thin fireshell in the fully radiative case (solid line) and in the adiabatic case (dotted line) 
starting from y = 10 2 and the same initial conditions as GRB 991216 [45], In the lower panel are plotted the corresponding values 
of the "effective" power-law index a e ff (see Eq.(l 1)), which is clearly not constant, is highly varying and systematically lower than 
the constant values 3 and 3/2 purported in the current literature (horizontal thin dotted lines). 

even At a ~ 10 3 s. For B < 3.0 x 10~ 4 , instead, At a presents a behavior which qualitatively follows the opposite of y 
(see middle panel of Fig. 4). 

Finally, in Fig. 8 we present three theoretical extended afterglow bolometric light curves together with the corre- 
sponding P-GRB peak luminosities for three different values of B. The duration of the P-GRBs has been assumed to 
be the same in the three cases (i.e. 5 s). The computations have been performed assuming the same E f °l and the same 
detailed CBM density profile of GRB 991216 [1]. In this picture we clearly see how, for B decreasing, the extended 
afterglow light curve "squeezes" itself on the P-GRB and the P-GRB peak luminosity increases. 

The "prompt emission" light curves of many GRBs present a small pulse preceding the main GRB event, with a 
lower peak flux and separated by this last one by a quiescent time. The nature of such GRB "precursors" is one of the 
most debated issues in the current literature [see e.g. Refs. 55, 32]. Already in 2001 [5], within the "fireshell" model, 
we proposed that GRB "precursors" are the P-GRBs emitted when the fireshell becomes transparent, and we gave 
precise estimates of the time sequence and intensities of the P-GRB and of the extended afterglow, recalled above. 

The radiation viewed in the co-moving frame of the accelerated baryonic matter is assumed to have a thermal 
spectrum and to be produced by the interaction of the CBM with the front of the expanding baryonic shell [49]. In 
Bernardini et al. [54] it is shown that, although the instantaneous spectrum in the co-moving frame of the optically thin 
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Figure 7. For 4 different values of E* e ±, we plot as a function of B the arrival time separation Af a between the P-GRB and the 
peak of the extended afterglow (i.e. the "quiescent time between the "precursor" and the main GRB event), measured in the source 
cosmological rest frame. This computation has been performed assuming a constant CBM density n chm = 1.0 particles/cm 3 . The 
points represents the actually numerically computed values, connected by straight line segments. 

fireshell is thermal, the shape of the final instantaneous spectrum in the laboratory frame is non-thermal. In fact, as 
explained in Ruffini et al. [49], the temperature of the fireshell is evolving with the co-moving time and, therefore, each 
single instantaneous spectrum is the result of an integration of hundreds of thermal spectra with different temperature 
over the corresponding EQTS. This calculation produces a non thermal instantaneous spectrum in the observer frame 
[54]. Another distinguishing feature of the GRBs spectra which is also present in these instantaneous spectra is the 
hard to soft transition during the evolution of the event [56, 43, 57, 58]. In fact the peak of the energy distributions 
E p drift monotonically to softer frequencies with time [54]. This feature explains the change in the power-law low 
energy spectral index [59] a which at the beginning of the prompt emission of the burst (t% = 2 s) is a = 0.75, and 
progressively decreases for later times [54]. In this way the link between E p and a identified in Crider et al. [56] is 
explicitly shown. 

The time-integrated observed GRB spectra show a clear power-law behavior. Within a different framework [see 
e.g. Ref. 60, and references therein] it has been argued that it is possible to obtain such power-law spectra from a 
convolution of many non power-law instantaneous spectra monotonically evolving in time. This result was recalled 
and applied to GRBs [61] assuming for the instantaneous spectra a thermal shape with a temperature changing with 
time. It was shown that the integration of such energy distributions over the observation time gives a typical power-law 
shape possibly consistent with GRB spectra. 

Our specific quantitative model is more complicated than the one considered by Blinnikov et al. [61]: the instan- 
taneous spectrum here is not a black body. Each instantaneous spectrum is obtained by an integration over the cor- 
responding EQTS: [47, 42] it is itself a convolution, weighted by appropriate Lorentz and Doppler factors, of ~ 10 6 
thermal spectra with variable temperature. Therefore, the time-integrated spectra are not plain convolutions of thermal 
spectra: they are convolutions of convolutions of thermal spectra [49, 54]. In Fig. 9 we present the photon number 
spectrum N(E) time-integrated over the 20 s of the whole duration of the prompt event of GRB 031203 observed 
by INTEGRAL [62]: in this way we obtain a typical non-thermal power-law spectrum which results to be in good 
agreement with the INTEGRAL data [62, 54] and gives a clear evidence of the possibility that the observed GRBs 
spectra are originated from a thermal emission [54] 
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Figure 8. We plot three theoretical extended afterglow bolometric light curves together with the corresponding P-GRB peak 
luminosities (the horizontal segments). The computations have been performed assuming the same E'± and CBM structure of GRB 
991216 and three different values of B. The P-GRBs have been assumed to have the same duration in the three cases, i.e. 5 s. For B 
decreasing, the extended afterglow light curve squeezes itself on the P-GRB. 



Before closing, we like to mention that, using the diagrams represented in Figs. 4-7, in principle one can compute 
the two free parameters of the fireshell model, namely E'J± and B, from the ratio between the total energies of the 
P-GRB and of the extended afterglow and from the temporal separation between the peaks of the corresponding 
bolometric light curves. None of these quantities depends on the cosmological model. Therefore, one can in principle 
use this method to compute the GRBs' intrinsic luminosity and make GRBs the best cosmological distance indicators 
available today. The increase of the number of observed sources, as well as the more accurate knowledge of their 
CBM density profiles, will possibly make viable this procedure to test cosmological parameters, in addition to the 
Amati relation [18, 36]. 



GRB 060607A: A COMPLETE ANALYSIS OF THE PROMPT EMISSION AND X-RAY 

FLARES. 

GRB 060607A is a very distant [z = 3.082, see Ref. 63] and energetic event [E iso ~ 10 53 erg, see Ref. 64]. Its BAT light 
curve shows a double-peaked structure with a duration of Tgo = (100 ±5) s [65]. The time-integrated spectrum over 
the T90 is best fit with a simple power-law model with an index T = 1.45 ±0.08 (Guidorzi, private communication). 
The XRT light curve shows a prominent flaring activity (at least three flares) superimposed to the normal afterglow 
decay [66]. 

The GRB 060607A main peculiarity is that the peak of the near-infrared (NIR) afterglow has been observed with the 
REM robotic telescope [64]. Interpreting the NIR light curve as corresponding to the afterglow onset as predicted by 
the fireball forward shock model [67, 28], it is possible to infer the initial Lorentz gamma factor of the emitting system 
that results to be F Q <~ 400 [64, 68, 69]. Moreover, these measurements seem to be consistent with an interstellar 
medium environment, ruling out the wind-like medium [64, 69]. 

We analyze GRB 060607A within the fireshell model [4, 5, 3]. We show that within this interpretation the N(E) 
spectrum of the prompt emission can be fitted in a satisfactory way by a convolution of thermal spectra as predicted 
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Figure 9. Three theoretically predicted time-integrated photon number spectra N(E), computed for GRB 031203 [54], are here 
represented for < t% < 5 s, 5 < t% < 10 s and 10 < t% < 20 s (dashed and dotted curves), where t% is the photon arrival time at 
the detector [4, 54]. The hard to soft behavior is confirmed. Moreover, the theoretically predicted time-integrated photon number 
spectrum N(E) corresponding to the first 20 s of the "prompt emission" (black bold curve) is compared with the data observed by 
INTEGRAL [62]. This curve is obtained as a convolution of 108 instantaneous spectra, which are enough to get a good agreement 
with the observed data. Details in Bernardini et al. [54]. 



by the model we applied [49, 50, 54]. The theoretical spectrum and light curve in the BAT energy band obtained are 
in good agreement with the observations, enforcing the plausibility of our approach. Moreover, in analogy with the 
case of GRB 011121 [70], we propose an interpretation of the observed X-ray flares as produced by overdense CBM 
clouds, in analogy with the gamma-ray light curve. 

In this preliminary analysis we deal only with the BAT and XRT observations, which are the basic contribution to 
the afterglow emission according to the fireshell model. We do not deal with the infrared emission that, on the contrary, 
is used in the current literature to estimate the dynamical quantities of the fireball in the forward external shock regime. 
Nevertheless, the initial value of Lorentz gamma factor we predict is compatible with the one deduced from the REM 
observations even under very different assumptions. 



GRB 060607A prompt emission 

Light curves 



In Fig. 10 we present the theoretical fit of Swift BAT light curves in different energy bands (15-25 keV, 25-50 keV, 
50-100 keV, 100-150 keV) of GRB 060607A. We identify the whole prompt emission with the peak of the extended 




Figure 10. Swift BAT (15-25 keV, 25-50 keV, 50-100 keV, 100-150 keV) light curves (points) compared with the theoretical 
ones (solid lines). 



afterglow emission, and the remaining part of the light curve with the decaying tail of the extended afterglow, according 
to our "canonical GRB" scenario [5, 3]. The temporal variability of the light curves has been reproduced assuming 
overdense spherical CBM regions [25]. The detailed structure of the CBM adopted is presented in Fig. 11. 

We therefore obtain for the two parameters characterizing the source in our model E'J± = 2.5 x 10 53 erg and 
B = 3.0 x 10~ 3 . This implies an initial e ± plasma with a total number of e^ pairs N e ± = 2.6 x 10 58 and an initial 
temperature T = 1.7 MeV. The theoretically estimated total isotropic energy emitted in the P-GRB is £p-G«B = 
1 .9%E'± = 4.7 x 10 51 erg, hence the P-GRB results to be undetectable if we assume a duration At p - gr i, > 10 s. 

After the transparency point at ro = 1.4 x 10 14 cm from the progenitor, the initial Lorentz gamma factor of the 
fireshell is y = 328. This value has been obtained adopting the exact solutions of the equations of motions of the 
fireshell [45] and using as initial conditions the two free parameters (E'J± and B) estimated from the simultaneous 
analysis of the BAT and XRT light curves. 



Time-integrated spectra 



We turn now to the analysis of the GRB 060607A prompt emission time-integrated spectrum. As discussed in 
previous works [49, 50, 54], even if the fireshell model assumes that the GRB spectrum is thermal in the comoving 
frame, the shape of the final spectrum in the laboratory frame is clearly non-thermal. In fact each single instantaneous 
spectrum is the result of a convolution of thermal spectra. In fact photons observed at the same arrival time are emitted 
at different comoving time [the so called EQTS, see Refs. 47, 42], hence with different temperatures. This calculation 
produces a non-thermal instantaneous spectrum in the observer frame. This effect is enhanced if we calculate the time- 
integrated spectrum: we perform two different integrations, one on the observation time and one on the EQTS, and 
what we get is a typical non-thermal power-law spectrum which results to be in good agreement with the observations 
(see Fig. 12). 
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Figure 11. Detailed structure of the CBM adopted: particle number density n c i, m (upper panel) and fraction effective emitting 
area 3% (lower panel) versus distance from the progenitor. The two X-ray flares corresponds in the upper panel to the huge increases 
in the CBM density that departs from the roughly power-law decrease observed. In the lower panel, the X-ray flares produce an 
increase of the emitting area which is not real but due to the lack of a complete 3-dimensional treatment of the interaction between 
the fireshell and the CBM (see text). 



The X-ray flares. 



We analyze now the X-ray flares observed by Swift XRT (0.2 — 10 keV) in the early part of the decaying phase 
of the extended afterglow. According to the fireshell model these flares have the same origin of the prompt emission, 
namely they are produced by the interaction of the fireshell with overdense CBM. As we can see in the upper panel 
in Fig. 13, the result obtained is compatible with the observations only for the second flare but not for the first one 
since its duration is longer. This discrepancy is due to the simple modeling adopted, namely the CBM is arranged in 
spherical shells [25]. This approximetion fails when the visible area of the fireshell is comparable with the size of the 
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Figure 12. Theoretically predicted time-integrated photon number spectrum N(E) corresponding to the — 15 s (upper panel), 
15 — 50 s (middle panel), and to the whole duration (Too = 100 s, lower panel) of the prompt emission (solid lines) compared with 
the observed spectra integrated in the same intervals. 
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Figure 13. Swift XRT (0.2-10 keV) light curve compared with, respectively, the theoretical one obtained assuming the CBM 
distributed in spherical shells (upper panel), the theoretical one obtained imposing a finite transverse dimension for the CBM cloud 
(middle panel) and the same theoretical curve in logarithmic scale (lower panel). 



CBM clouds. 

To solve this problem, following the results obtained for GRGB 011121 [70], we tried to account for the three- 
dimensional structure of the CBM clouds by "cutting" the emission at a certain angle Q c ioud from the line of sight, 
corresponding to the transverse dimension of the CBM cloud, until the duration of the flare 8t / t tot is compatible with 
the observation (see Fig. 13 middle and lower panels). It is worth to observe that with this procedure we kee the value 
of R constant during the flare. Hence the increase in R that we obtained in our previous analysis (see Fig. 1 1) is not 
real but it compensates the fact that spherical approximation is not valid at this stage. 

This procedure affects the dynamics of the fireshell, so the light curve after the "cut" is meaningless. Nevertheless, 
it is a confirmation that it is possible to obtain arbitrarily short flares by the interaction with the CBM. 



Conclusions 

We presented the analysis of GRB 060607A within the fireshell model [4, 5, 3]. According to the "canonical GRB" 
scenario [5, 3] we interpreted the whole prompt emission as the peak of the extended afterglow emission, and the 
remaining part of the light curve with the decaying tail of the extended afterglow. We found in this second case that 
the P-GRB is too faint to be detected, as we expected from our interpretation. The theoretical light curves obtained are 
well in agreement with the observations in all the Swift BAT energy bands. 

Furthermore, the initial Lorentz gamma factor of the fireshell, obtained adopting the exact solutions of its equations 
of motions [45] and as initial condition the free parameters of the fireshell estimated by the simultaneous analysis 
of the BAT and XRT light curves, is y = 328. In this preliminary analysis we deal only with the BAT and XRT 
observations, which are the basic contribution to the afterglow emission according to the fireshell model. We do not 
deal with the infrared emission that, on the contrary, is used in the current literature to estimate the dynamical quantities 
of the fireball in the forward external shock regime. Nevertheless, the initial value of Lorentz gamma factor we predict 
is compatible with the one deduced from the REM observations even under very different assumptions, r o <~ 400 
[64,68,69]. 

We investigated also the GRB 060607A prompt emission spectra integrated in different time intervals assuming 
a thermal spectrum in the comoving frame. The results obtained show clearly that, after the correct space-time 
transformations, both the instantaneous and the time-integrated spectra in the observer frame have nothing to do with 
a Planckian distribution, but they have a power-law shape, thus confirming our previous analyses [54, 71]. 

Finally we analyzed the X-ray flares observed by Swift XRT (0.2 — 10 keV) in the early part of the decaying phase 
of the extended afterglow. According to the fireshell model these flares have the same origin of the prompt emission, 
namely they are produced by the interaction of the fireshell with overdense CBM. We found that our theoretical light 
curve is not compatible with the observations since in such regime the one-dimensional approximation fails. Following 
the results obtained for GRGB 011121 [70], we tried to account for the three-dimensional structure of the CBM clouds 
by "cutting" the emission at a certain angle 9 c i„ uc i from the line of sight, corresponding to the transverse dimension of 
the CBM cloud. We obtain in this way a flare whose duration 8t /t tot is compatible with the observation. 

THE NORRIS & BONNEL KIND OF SOURCES: THE NEW CLASS OF "FAKE - 

DISGUISED" SHORT GRBS 

We now present the theoretical understanding, within the fireshell model, of a new class of sources, pioneered by 
Norris and Bonnell [7]. This class is characterized by an occasional softer extended emission after an initial spikelike 
emission. The softer extended emission has a peak luminosity smaller than the one of the initial spikelike emission. 
This has misled the understanding of the correct role of the extended afterglow. As shown in the prototypical case 
of GRB 970228 [see below and Ref. 33], the initial spikelike emission can be identified with the P-GRB and the 
softer extended emission with the peak of the extended afterglow. Crucial is the fact that the time-integrated extended 
afterglow luminosity is much larger than the P-GRB one, and this fact unquestionably identifies GRB 970228 as a 
canonical GRB with B > 10~ 4 . The consistent application of the fireshell model allowed to compute the CBM porosity, 
filamentary structure and average density which, in that specific case, resulted to be n c b m ~ 10~ 3 particles/cm 3 [33]. 
This explained the peculiarity of the low extended afterglow peak luminosity and of its much longer time evolution. 
These features are not intrinsic to the progenitor nor to the black hole, but they uniquely depend on the peculiarly 
low value of the CBM density, typical of galactic halos. If one takes the same total energy, baryon loading and CBM 
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Figure 14. A sketch summarizing the "canonical GRB" scenario. 



distribution as in GRB 970228, and rescales the CBM density profile by a constant numerical factor in order to raise its 
average value from 10~ 3 to 1 particles/cm 3 , he obtains a GRB with a much larger extended afterglow peak luminosity 
and a much reduced time scale. Such a GRB would appear a perfect traditional "long" GRB following the current 
literature [see below and Ref. 33]. This has led us to expand the traditional classification of GRBs to three classes: 
"genuine" short GRBs, "fake" or "disguised" short GRBs, and all the remaining "canonical" ones [see Fig. 14 and 
Ref. 72]. 



GRB 970228 and a class of GRBs with an initial spikelike emission 

GRB 970228 was detected by the Gamma-Ray Burst Monitor (GRBM, 40-700 keV) and Wide Field Cameras 
(WFC, 2-26 keV) on board BeppoSAX on February 28.123620 UT [73]. The burst prompt emission is characterized 
by an initial 5 s strong pulse followed, after 30 s, by a set of three additional pulses of decreasing intensity [73]. Eight 
hours after the initial detection, the NFIs on board BeppoSAX were pointed at the burst location for a first target 
of opportunity observation and a new X-ray source was detected in the GRB error box: this is the first "afterglow" 
ever detected [74]. A fading optical transient has been identified in a position consistent with the X-ray transient [75], 
coincident with a faint galaxy with redshift z — 0.695 [76] . Further observations by the Hubble Space Telescope clearly 
showed that the optical counterpart was located in the outskirts of a late-type galaxy with an irregular morphology [77] . 

The BeppoSAX observations of GRB 970228 prompt emission revealed a discontinuity in the spectral index 
between the end of the first pulse and the beginning of the three additional ones [74, 73, 57]. The spectrum during the 
first 3 s of the second pulse is significantly harder than during the last part of the first pulse [73, 57], while the spectrum 
of the last three pulses appear to be consistent with the late X-ray afterglow [73, 57]. This was soon recognized by 
Frontera et al. [73, 57] as pointing to an emission mechanism producing the X-ray afterglow already taking place after 
the first pulse. 

The simultaneous occurrence of an extended afterglow with total time-integrated luminosity larger than the P-GRB 
one, but with a smaller peak luminosity, is indeed explainable in terms of a peculiarly small average value of the CBM 
density and not due to the intrinsic nature of the source. In this sense, GRBs belonging to this class are only "fake" or 
"disguised" short GRBs. We show that GRB 970228 is a very clear example of this situation. We identify the initial 
spikelike emission with the P-GRB, and the late soft bump with the peak of the extended afterglow. GRB 970228 
shares the same morphology and observational features with the sources analyzed by Norris and Bonnell [7] as well as 
with e.g. GRB 050709 [78], GRB 050724 [79] and GRB 060614 [see next section and Ref. 80]. Therefore, we propose 
GRB 970228 as a prototype for this new GRB class. 
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Figure 15. The "canonical GRB" light curve theoretically computed for the prompt emission of GRB 970228. BeppoSAX GRBM 
(40-700 keV, above) and WFC (2-26 keV, below) light curves (data points) are compared with the extended afterglow peak 
theoretical ones (solid lines). The onset of the extended afterglow coincides with the end of the P-GRB (represented qualitatively 
by the dotted lines). For this source we have B ~ 5.0 x 10~ 3 and (n c j m ) ~ 10~ 3 particles/cm 3 . Details in Bernardini et al. [33, 81]. 



The analysis of GRB 970228 prompt emission 



In Fig. 15 we present the theoretical fit of BeppoSAX GRBM (40-700 keV) and WFC (2-26 keV) light curves of 
GRB 970228 prompt emission [73]. Within our "canonical GRB" scenario we identify the first main pulse with the P- 
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Figure 16. The CBM density profile we assumed to reproduce the last three pulses of the GRB 970228 prompt emission (red 
line), together with its average value (n c i, m ) = 9.5 x 10~ 4 particles/cm 3 (green line). 



GRB and the three additional pulses with the extended afterglow peak emission, consistently with the above mentioned 
observations by Costa et al. [74] and Frontera et al. [73]. Such last three pulses have been reproduced assuming three 
overdense spherical CBM regions (see Fig. 16) with a very good agreement (see Fig. 15). 

We therefore obtain for the two parameters characterizing the source in our model E'f = 1.45 x 10 54 erg and 
B = 5.0 x 10~ 3 . This implies an initial e ± plasma created between the radii r\ = 3.52 x 10 7 cm and r2 = 4.87 x 10 8 cm 
with a total number of e ± pairs N e ± = 1 .6 x 10 59 and an initial temperature T = 1 .7 MeV. The theoretically estimated 
total isotropic energy emitted in the P-GRB is £/>_gj?b = \.\%E to l = 1.54 x 10 52 erg, in excellent agreement with the 
one observed in the first main pulse (Ep^ GRB <~ 1 .5 x 10 52 erg in 2 — 700 keV energy band, see Fig. 15), as expected 
due to their identification. After the transparency point at ro = 4.37 x 10 14 cm from the progenitor, the initial Lorentz 
gamma factor of the fireshell is yo = 199. On average, during the extended afterglow peak emission phase we have for 
the CBM {M) = 1 .5 x 10~ 7 and (n cbm ) = 9.5 x 10~ 4 particles/cm 3 . This very low average value for the CBM density 
is compatible with the observed occurrence of GRB 970228 in its host galaxy's halo [77, 75, 82] and it is crucial in 
explaining the light curve behavior. 

The values of E'J± and B we determined are univocally fixed by two tight constraints. The first one is the total energy 
emitted by the source all the way up to the latest extended afterglow phases (i.e. up to <~ 10 6 s). The second one is 
the ratio between the total time-integrated luminosity of the P-GRB and the corresponding one of the whole extended 
afterglow (i.e. up to <~ 10 6 s). In particular, in GRB 970228 such a ratio results to be <~ 1.1% (see Fig. 5). However, 
the P-GRB peak luminosity actually results to be much more intense than the extended afterglow one (see Fig. 15). 
This is due to the very low average value of the CBM density (n c i, m ) = 9.5 x 10~ 4 particles/cm 3 , which produces a 
less intense extended afterglow emission. Since the extended afterglow total time-integrated luminosity is fixed, such 
a less intense emission lasts longer than what we would expect for an average density (n c i, m ) <~ 1 particles/cm 3 . 



Rescaling the CBM density 



We present now an explicit example in order to probe the crucial role of the average CBM density in explaining 
the relative intensities of the P-GRB and of the extended afterglow peak in GRB 970228. We keep fixed the basic 
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Figure 17. The theoretical fit of the BeppoSAX GRBM observations (solid line, see Fig. 15) is compared with the extended 
afterglow light curve in the 40-700 keV energy band obtained rescaling the CBM density to (n c i, m ) = 1 particle/cm 3 keeping 
constant its shape and the values of the fundamental parameters of the theory E x °£ and B (double dotted line). The P-GRB duration 
and luminosity (dotted line), depending only on E*°l and B, are not affected by this process of rescaling the CBM density. 



parameters of the source, namely the total energy E'± and the baryon loading B, therefore keeping fixed the P-GRB 
and the extended afterglow total time-integrated luminosities. Then we rescale the CBM density profile given in Fig. 
16 by a constant numerical factor in order to raise its average value to the standard one iiiisni) — 1 particle/cm 3 . We 
then compute the corresponding light curve, shown in Fig. 17. 

We notice a clear enhancement of the extended afterglow peak luminosity with respect to the P-GRB one in 
comparison with the fit of the observational data presented in Fig. 15. The two light curves actually crosses at 

~ 1.8 x 10 4 s since their total time-integrated luminosities must be the same. The GRB "rescaled" to (n, jm ) = 1 
particle/cm 3 appears to be totally similar to, e.g., GRB 050315 [71] and GRB 991216 [1, 49, 2]. 

It is appropriate to emphasize that, although the two underlying CBM density profiles differ by a constant numerical 
factor, the two extended afterglow light curves in Fig. 17 do not. This is because the absolute value of the CBM 
density at each point affects in a non-linear way all the following evolution of the fireshell due to the feedback on 
its dynamics [45]. Moreover, the shape of the surfaces of equal arrival time of the photons at the detector (EQTS) 
is strongly elongated along the line of sight [42]. Therefore photons coming from the same CBM density region are 
observed over a very long arrival time interval. 



GRB 970228 and the Amati relation 

We turn now to the "Amati relation" [13, 20] between the isotropic equivalent energy emitted in the prompt emission 
Ei S0 and the peak energy of the corresponding time-integrated spectrum E p j in the source rest frame. It has been shown 
by Amati et al. [13], Amati [20] that this correlation holds for almost all the "long" GRBs which have a redshift and 
an E p j measured, but not for the ones classified as "short" [20]. If we focus on the "fake" or "disguised" short GRBs, 
namely the GRBs belonging to this new class, at least in one case [GRB 050724, see Ref. 79] it has been shown that 
the correlation is recovered if also the extended emission is considered [83]. 

It clearly follows from our treatment that for the "canonical GRBs" with large values of the baryon loading and high 
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Figure 18. The estimated values for E p j and Ej so obtained by our analysis (black dot) compared with the "Amati relation" [13]: 
the solid line is the best fitting power law [20] and the dashed lines delimit the region corresponding to a vertical logarithmic 
deviation of 0.4 [20]. The uncertainty in the theoretical estimated value for E p j has been assumed conservatively as 20%. 



{ n cbm)i which presumably are most of the GRBs for which the correlation holds, the leading contribution to the prompt 
emission is the extended afterglow peak emission. The case of the "fake" or "disguised" short GRBs is completely 
different: it is crucial to consider separately the two components since the P-GRB contribution to the prompt emission 
in this case is significant. 

To test this scenario, we evaluated from our fit of GRB 970228 Ej S0 and E p j only for the extended afterglow peak 
emission component, i.e. from = 37 s to t% = 81.6 s. We found an isotropic energy emitted in the 2-400 keV 
energy band Ej S0 = 1 .5 x 10 52 erg, and E p j = 90.3 keV. As it is clearly shown in Fig. 18, the sole extended afterglow 
component of GRB 970228 prompt emission is in perfect agreement with the Amati relation. If this behavior is 
confirmed for other GRBs belonging to this new class, this will enforce our identification of the "fake" or "disguised" 
short GRBs. This result will also provide a theoretical explanation for the the apparent absence of such correlation for 
the initial spikelike component in the different nature of the P-GRB. 



Conclusions 

We conclude that GRB 970228 is a "canonical GRB" with a large value of the baryon loading quite near to the 
maximum B ~ 10~ 2 (see Fig. 5). The difference with e.g. GRB 050315 [71] or GRB 991216 [1, 49, 2] is the low 
average value of the CBM density (n c & m ) <~ 10~ 3 particles/cm 3 which deflates the extended afterglow peak luminosity. 
Hence, the predominance of the P-GRB, coincident with the initial spikelike emission, over the extended afterglow 
is just apparent: 98.9% of the total time-integrated luminosity is indeed in the extended afterglow component. Such 



a low average CBM density is consistent with the occurrence of GRB 970228 in the galactic halo of its host galaxy 
[77, 75], where lower CBM densities have to be expected [82]. 

We propose GRB 970228 as the prototype for the new class of GRBs comprising GRB 060614 and the GRBs 
analyzed by Norris and Bonnell [7]. We naturally explain the hardness and the absence of spectral lag in the initial 
spikelike emission with the physics of the P-GRB originating from the gravitational collapse leading to the black 
hole formation. The hard-to-soft behavior in the extended afterglow is also naturally explained by the physics of the 
relativistic fireshell interacting with the CBM, clearly evidenced in GRB 031203 [54] and in GRB 050315 [71]. Also 
justified is the applicability of the Amati relation to the sole extended afterglow component [see Refs. 20, 83]. 

This class of GRBs with z ~ 0.4 appears to be nearer than the other GRBs detected by Swift [z ~ 2.3, see Ref. 
84]. This may be explained by the extended afterglow peak luminosity deflation. The absence of a jet break in those 
afterglows has been pointed out [79, 85], consistently with our spherically symmetric approach. Their association with 
non-star-forming host galaxies appears to be consistent with the merging of a compact object binary [86, 87]. It is here 
appropriate, however, to caution on this conclusion, since the association of GRB 060614 and GRB 970228 with the 
explosion of massive stars is not excluded [88, 89]. 

Most of the sources of this class appear indeed not to be related to bright "Hypernovae", to be in the outskirts of their 
host galaxies [87, see above] and a consistent fraction of them are in galaxy clusters with CBM densities (n c f, m ) <~ 10~ 3 
particles/cm 3 [see e.g. Ref. 90, 91]. This suggests a spiraling out binary nature of their progenitor systems [92] made 
of neutron stars and/or white dwarfs leading to a black hole formation. 

Moreover, we verified the applicability of the Amati relation to the sole extended afterglow component in GRB 
970228 prompt emission, in analogy with what happens for some of the GRBs belonging to this new class. In fact it 
has been shown by Amati [20, 83] that the "fake" or "disguised" short GRBs do not fulfill the E p j-Ei S0 correlation 
when the sole spiklike emission is considered, while they do if the long soft bump is included. Since the spikelike 
emission and the soft bump contributions are comparable, it is natural to expect that the soft bump alone will fulfill the 
correlation as well. 

Within our "canonical GRB" scenario the sharp distinction between the P-GRB and the extended afterglow provide 
a natural explanation for the observational features of the two contributions. We naturally explain the hardness and the 
absence of spectral lag in the initial spikelike emission with the physics of the P-GRB originating from the gravitational 
collapse leading to the black hole formation. The hard-to-soft behavior in the extended afterglow is also naturally 
explained by the physics of the relativistic fireshell interacting with the CBM, clearly evidenced in GRB 031203 [54] 
and in GRB 050315 [71]. Therefore, we expect naturally that the E p j-Ei SO correlation holds only for the extended 
afterglow component and not for the P-GRB. Actually we find that the correlation is recovered for the extended 
afterglow peak emission of GRB 970228. 

In the original work by Amati et al. [13], Amati [20] only the prompt emission is considered and not the late 
afterglow one. In our theoretical approach the extended afterglow peak emission contributes to the prompt emission 
and continues up to the latest GRB emission. Hence, the meaningful procedure within our model to recover the Amati 
relation is to look at a correlation between the total isotropic energy and the peak of the time-integrated spectrum of the 
whole extended afterglow. A first attempt to obtain such a correlation has already been performed using GRB 050315 
as a template, giving very satisfactory results (see section Theoretical background for GRBs' empirical correlations). 

THE "FIRESHELL" MODEL AND GRB PROGENITORS 

"Long" GRBs are traditionally related in the current literature to the idea of a single progenitor, identified as a 
"collapsar" [93]. Similarly, short GRBs are assumed to originate from binary mergers formed by white dwarfs, neutron 
stars, and black holes in all possible combinations [see e.g. Refs. 94, 95, 96, 97, 27, 28, and references therein]. It 
has been also suggested that short and long GRBs originate from different galaxy types. In particular, short GRBs 
are proposed to be associated with galaxies with low specific star forming rate [see e.g. Ref. 98]. Some evidences 
against such a scenario have been however advanced, due to the small sample size and the different estimates of the 
star forming rates [see e.g. Ref. 99]. However, the understanding of GRB structure and of its relation to the CBM 
distribution, within the fireshell model, leads to a more complex and interesting perspective than the one in the current 
literature. 

The first general conclusion of the "fireshell" model [5] is that, while the time scale of "short" GRBs is indeed 
intrinsic to the source, this does not happen for the "long" GRBs: their time scale is clearly only a function of the 
instrumental noise threshold. This has been dramatically confirmed by the observations of the Swift satellite [see Fig. 
19 and Ref. 101]. Among the traditional classification of "long" GRBs we distinguish two different sub-classes of 




Figure 19. The theoretical light curves in the 15 — 150 keV (solid line) and 0.2— 10 kev (dotted line) energy bands compared 
with XRT observations of GRB 050315 [100]. The horizontal dashed lines correspond to different possible instrumental thresholds. 
It is clear that long GRB durations are just functions of the observational threshold. Details in Ruffini et al. [101]. 



events, none of which originates from collapsars. 

The first sub-class contains "long" GRBs particularly weak (Ej S0 <~ 10 50 erg) and associated with Supernovae (SNe) 
Ib/c. In fact, it has been often proposed that such GRBs, only observed at smaller redshift 0.0085 < z < 0.168, form 
a different class, less luminous and possibly much more numerous than the high luminosity GRBs at higher redshift 
[102, 103, 104, 105]. Therefore in the current literature they have been proposed to originate from a separate class 
of progenitors [106, 107]. Within our "fireshell" model, they originate in a binary system formed by a neutron star, 
close to its critical mass, and a companion star, evolved out of the main sequence. They produce GRBs associated with 
SNe Ib/c, via the "induced gravitational collapse" process [6]. The low luminosity of these sources is explained by the 
formation of a black hole with the smallest possible mass: the one formed by the collapse of a just overcritical neutron 
star [108, 51]. 

A second sub-class of "long" GRBs originates from merging binary systems, formed either by two neutron stars or 
a neutron star and a white dwarf. A prototypical example of such systems is GRB 970228. The binary nature of the 
source is inferred by its migration from its birth location in a star forming region to a low density region within the 
galactic halo, where the final merging occurs [33]. The location of such a merging event in the galactic halo is indeed 
confirmed by optical observations of the GRB 970228 afterglow [77, 75]. The crucial point is that, as recalled above, 
GRB 970228 is a "canonical" GRB with B > 10~ 4 "disguised" as a short GRB. We are going to see in the following 
that GRB 060614 also comes from such a progenitor class. 

If the binary merging would occur in a region close to its birth place, with an average density of 1 particle/cm 3 , 
the GRB would appear as a traditional high-luminosity "long" GRB, of the kind currently observed at higher redshifts 
(see above, Fig. 17), similar to, e.g., GRB 050315 [71]. 



Within our approach, therefore, there is the distinct possibility that all GRB progenitors are formed by binary 
systems, composed by neutron stars, white dwarfs, or stars evolved out of the main sequence, in different combinations. 
The case of the "genuine" short GRBs is currently being examined within the "fireshell" model. 



GRB 060614: a "fake" or "disguised" short GRB from a merging binary system 

GRB 060614 [80, 109] has imposed to the general attention of the Gamma-Ray Burst's (GRB's) scientific commu- 
nity because it is the first clear example of a nearby (z = 0. 125), long GRB not associated with a bright Ib/c Supernova 
(SN) [88, 110]. It has been estimated that, if present, the SN-component should be about 200 times fainter than the 
archetypal SN 1998bw associated to GRB 980425; moreover, it would also be fainter (at least 30 times) than any 
stripped-envelope SN ever observed [111]. 

Within the standard scenario, long duration GRBs (T90 > 2 s) are thought to be produced by SN events during the 
collapse of massive stars in star forming regions ["collapsar", see Ref. 93]. The observations of broad-lined and bright 
type Ib/c SNe associated with GRBs are often reported to favor this scenario [see Ref. 112, and references therein]. 
The ansatz has been advanced that every long GRB should have a SN associated with it [113]. Consequently, in all 
nearby long GRBs (z < 1) the SN emission should be observed. 

For these reasons the case of GRB 060614 is indeed revolutionary. Some obvious hypothesis have been proposed 
and ruled out: the chance superposition with a galaxy at low redshift [110] and the strong dust obscuration and 
extinction [1 14]. Appeal has been made to the possible occurrence of an unusually low luminosity stripped-envelope 
core-collapse SN [88]. 

The second novelty of GRB 060614 is that it challenges the traditional separation between Long Soft GRBs and 
Short Hard GRBs. Traditionally [8, 9], the "short" GRBs have Tgo < 2 s, present an harder spectrum and negligible 
spectral lag, and are assumed to originate from merging of two compact objects, i.e. two neutron stars or a neutron 
star and a black hole [see e.g. Ref. 94, 95, 96, 97, 27, 28, and references therein]. GRB 060614 lasts about one 
hundred seconds [?9o = (102 ±5) s; see Ref. 80], it fulfills the Ep St — Ei so correlation [115], and therefore it should 
be traditionally classified as a "long" GRB. However, its morphology is different from typical long GRBs, similar 
to the one of GRB 050724, traditionally classified as a short GRB [113, 116]. Its optical afterglow luminosity is 
intermediate between the traditional long and short ones [117]. Its host galaxy has a moderate specific star formation 
rate [Ruost ~ 2M s y~ 1 (L* ) ~ 1 , M V H„st ~ —15.5; see Refs. 114, 88]. The spectral lag in its light curves is very small or 
absent [80]. All these features are typical of the short GRBs. 

A third peculiarity of GRB 060614 is that its 15-150 keV light curve presents a short, hard and multi -peaked episode 
(about 5 s). Such an episode is followed by a softer, prolonged emission that manifests a strong hard to soft evolution 
in the first 400 s of data [109]. The total fluence in the 15-150 keV energy band is F = (2.17 ±0.04) x 10~ 5 erg/cm 2 , 
the 20% emitted during the initial spikelike emission, where the peak luminosity reaches the value of 300 keV before 
decreasing until 8 keV during the BAT-XRT overlap time (about 80 s). 

These apparent contradictions find a natural explanation in the framework of the "fireshell" model. Within the 
fireshell model, the occurrence of a GRB-SN is not a necessity. The origin of all GRBs is traced back to the formation 
of a black hole, either occurring in a single process of gravitational collapse, or in a binary system composed by a 
neutron star and a companion star evolved out of the main sequence, or in a merging binary system composed by 
neutron stars and/or white dwarfs in all possible combinations. The occurrence of a GRB-SN is indeed only one of the 
possibilities, linked, for example, to the process of "induced gravitational collapse" [6, 108, 51]. 

We here show how the "fireshell" model can explain all the above mentioned GRB 060614 peculiarities and solve 
the apparent contradictions. In doing so, we also infer constraints on the astrophysical nature of the GRB 060614 
progenitors. In turn, these conclusions lead to a new scenario for all GRBs. We can confirm a classification of GRBs 
in "genuine" short, "fake" or "disguised" short, and, finally, all the remaining "canonical" GRBs. The connection 
between this new classification and the nature of GRB progenitors is quite different from the traditional one in the 
current literature. 



The fit of the observed luminosity 

In this scenario, GRB 060614 is naturally interpreted as a "disguised" short GRB. We have performed the analysis 
of the observed light curves in the 15-150 keV energy band, corresponding to the y-ray emission observed by the BAT 



instrument on the Swift satellite, and in the 0.2-10 keV energy band, corresponding to the X-ray component from the 
XRT instrument on Swift satellite. We do not address in this paper the issue of the optical emission, that represent less 
than 10% of the total energy of the GRB. From this fit (see Figs. 20, 22) we have derived the total initial energy Ef ot , 
the value of B as well as the effective CBM distribution (see Fig. 21). We find Ef ot = 2.94 x 10 erg, that accounts for 
the bolometric emission of both the P-GRB and the extended afterglow. Such a value is compatible with the observed 
E iso ~ 2.5 x 10 51 erg [80]. The value of B is B = 2.8 x 10~ 3 , that corresponds to the lowest one of all the GRBs we have 
examined (see Fig. 5). It corresponds to a canonical GRB with a very clear extended afterglow predominance over the 
P-GRB. From the model, having determined Ef ot and B, we can compute the theoretical expected P-GRB energetics 
Ep-grb [5]. We obtain E P - GRB ~ 1.15 x 10 50 erg, that is in good agreement with the observed one E iso ^ p ~ 1.18 x 10 50 
erg [80]. The Lorentz Gamma Factor at the transparency results to be y = 346, one of the highest of all the GRBs we 
have examined. 

In Fig. 20 we plot the comparison between the BAT observational data of the GRB 0606014 prompt emission 
in the 15-150 keV energy range and the P-GRB and extended afterglow light curves computed within our model. 
The temporal variability of the extended afterglow peak emission is due to the inhomogeneities in the effective CBM 
density (see Figs. 20, 21). Toward the end of the BAT light curve, the good agreement between the observations and the 
fit is affected by the Lorentz gamma factor decrease and the corresponding increase of the maximum viewing angle. 
The source visible area becomes larger than the typical size of the filaments. This invalidates the radial approximation 
we use for the CBM description. To overcome this problem it is necessary to introduce a more detailed three- 
dimensional CBM description, in order to avoid an over-estimated area of emission and, correspondingly, to describe 
the sharpness of some observed light curves. We are still working on this issue [25, 118, 70, 1 19]. 

We turn now to the crucial determination of the CBM density, which is derived from the fit. At the transparency 
point it resulted to be n c b m = 4.8 x 10~ 3 particles/cm 3 (see Fig. 21). This density is compatible with the typical values 
of the galactic halos. During the peak of the extended afterglow emission the effective average CBM density decreases 
reaching (n c & m ) = 2.25 x 10~ 5 particles/cm 3 , possibly due to an occurring fragmentation of the shell [51] or due to a 
fractal structure in the CBM. The value resulted to be on average (&) = 1.72 x 10~ 8 . It is interesting to emphasize 
the striking analogy of the numerical value and the overall radial dependence of the CBM density in the present case 
of GRB 060614 when compared and contrasted with the ones of GRB 970228 [33]. 

Concerning the 0.2-10 keV light curve of the decaying phase of the afterglow, observed by the XRT instrument, we 
have also reproduced very satisfactorily both the hard decrease in the slope and the plateau of the light curve keeping 
constant the effective CBM density and changing only The result of this analysis is reported in Fig. 22. We assume 
in this phase n c f, m = 4.70 x 10~ 6 particles/cm -3 . The average value of the & parameter is {&) = 1.27 x 10~ 2 . The 
drastic enhancement in the & parameter with respect to the values at the peak of the extended afterglow is consistent 
with similar features encountered in other sources we have studied: GRB 060218 presents a bump of five orders of 
magnitude [51], in GRB 060710 the bump is of about four orders of magnitude (see Izzo et al., in preparation) while 
in GRB 050315 there is a three orders of magnitude bump [71]. In these last two cases, we find the occurrence of 
the enhancement of & between r=2 x 10 17 cm and r=3 x 10 17 cm, just like for GRB 060614, for which we have the 
bump at r=3.5 x 10 17 cm. The time of the bump approximately corresponds to the appearance of the optical emission 
observed in GRB 060614 and, more in general, to the onset of the second component of the Willingale et al. [120] 
scheme for GRBs. 



Conclusions 

GRB 060614 presents three major novelties, which challenges the most widespread theoretical models and which 
are strongly debated in the current literature. The first one is that it challenges the traditional separation between Long 
Soft GRBs and Short Hard GRBs [80]. The second one is that it presents a short, hard and multi-peaked episode, 
followed by a softer, prolonged emission with a strong hard to soft evolution [80, 109]. The third one is that it is the 
first clear example of a nearby, long GRB not associated with a bright SN Ib/c [88, 110]. All these three issues are 
naturally explained within our "fireshell" model, which allows a detailed analysis of the temporal behavior of the signal 
originating up to a distance r ~ 10 17 -10 18 cm from the black hole, and relates, with all the relativistic transformations, 
the arrival time to the CBM structure and the relativistic parameters of the fireshell. 

One of the major outcome of the Swift observation of, e.g., GRB 050315 [100, 71] has been the confirmation that 
long GRB duration is not intrinsic to the source but it is merely a function of the instrumental noise threshold [101]. 
GRB 060614 represents an additional fundamental progress in clarifying the role of the CBM density in determining 



8.0x10 



7.0x10' 




-1.0x10 



40 60 
Detector arrival time (s) 



120 



Figure 20. The BAT 15-150 keV light curve (points) at 1 s time resolution is compared with the corresponding theoretical 
extended afterglow light curve we compute (solid line). The onset of the extended afterglow is at the end of the P-GRB (qualitatively 
sketched in dashed lines and delimited by dashed tick vertical lines). Therefore the zero of the temporal axis is shifted by 5.5 s with 
respect to the BAT trigger time. The peaks of the extended afterglow light curves are labeled to match them with the corresponding 
CBM density peak in Fig. 21. In the upper right corner there is an enlargement of the P-GRB at 50ms time resolution [reproduced 
from Ref. 109], showing its structure. 



the GRB morphology. It confirms the results presented in GRB 970228 [33], that is the prototype of the new class 
of "fake" or "disguised" short GRBs. They correspond to canonical GRBs with an extended afterglow emission 
energetically predominant with respect to the P-GRB one and a baryon loading B > 10~ 4 . The sharp spiky emission 
corresponds to the P-GRB. As recalled above, a comparison of the luminosities of the P-GRB and of the extended 
afterglow is indeed misleading: it follows from the low average CBM density inferred from the fit of the fireshell 
model, which leads to n c f, m ~ 10~ 3 particles/cm 3 . Therefore such a feature is neither intrinsic to the progenitor nor to 
the black hole, but it is only indicative of the CBM density at the location where the final merging occurs. GRB 060614 
is a canonical GRB and it is what would be traditionally called a "long" GRB if it had not exploded in a specially low 
CBM density environment. GRB 060614 must necessarily fulfill, and indeed it does, the Amati relation. This happens 
even taking into account the entire prompt emission mixing together the P-GRB and the extended afterglow [115], due 
to the above recalled energetic predominance of the extended afterglow [see also Ref. 36]. These results justify the 
occurrence of the above mentioned first two novelties. 

The low value of the CBM density is compatible with a galactic halo environment. This result points to an old 
binary system as the progenitor of GRB 060614 and it justifies the above mentioned third novelty: the absence of 
an associated SN Ib/c [see also Ref. 121]. Such a binary system departed from its original location in a star forming 
region and spiraled out in a low density region of the galactic halo [see e.g. Ref. 92]. The energetic of this GRB is 
about two orders of magnitude smaller than the one of GRB 970228 [33]. A natural possible explanation is that instead 
of a neutron star - neutron star merging binary system we are in presence of a white dwarf - neutron star binary. We 
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Figure 21. Here are the plot of the effective CBM density (solid line) and of the 3% parameter (dotted line) versus the radial 
coordinate of the shell. The CBM density peaks are labeled to match them with the corresponding extended afterglow light curve 
peaks in Fig. 20. They corresponds to filaments of characteristic size Ar ~ 10 15 cm and density contrast An cbm /(n cbm ) ~ 20 
particles/cm 3 . 



therefore agree, for different reasons, with the identification proposed by Davies et al. [121] for the GRB 060614 
progenitor. In principle, the nature of the white dwarf, with typical radius on the order of 10 3 km, as opposed to the 
one of the neutron star, typically on the order of 10 km, may manifest itself in characteristic signatures in the structure 
of the P-GRB (see Fig. 20). 

It is interesting that these results lead also to three major new possibilities: 

• The majority of GRBs declared as shorts [see e.g. Ref. 116] are likely "disguised" short GRBs, in which the 
extended afterglow is below the instrumental threshold. 

• The observations of GRB 060614 offer the opportunity, for the first time, to analyze in detail the structure of a 
P-GRB lasting 5 s. This feature is directly linked to the physics of the gravitational collapse originating the GRB. 
Recently, there has been a crucial theoretical physics result, showing that the characteristic time constant for the 
thermalization for an plasma is on the order of 10~ 13 s [122]. Such a time scale still applies for an e ± plasma 
with a baryon loading on the order of the one observed in GRBs [40]. The shortness of such a time scale, as well 
as the knowledge of the dynamical equations of the optically thick phase preceding the P-GRB emission [31], 
implies that the structure of the P-GRB is a faithful representation of the gravitational collapse process leading 
to the formation of the black hole [123]. In this respect, it is indeed crucial that the Swift data on the P-GRB 
observed in GRB 060614 (see Fig. 20) appear to be highly structured all the way to time scale of 0.1 s. This 
opens a new field of research: the study of the P-GRB structure in relation to the process of gravitational collapse 
leading to the GRB. 

• If indeed the binary nature of the progenitor system and the peculiarly low CBM density n c & m ~ 10~ 3 
particles/cm 3 will be confirmed for all "fake" or "disguised" GRBs, then it is very likely that the traditionally 
"long" high luminosity GRBs at higher redshift also originates from the merging of binary systems formed by 
neutron stars and/or white dwarfs occurring close to their birth location in star forming regions with n c b m <~ 1 
particle/cm 3 (see Fig. 17). 




Detector arrival time (s) 

Figure 22. The XRT 0.2-10 keV light curve (points) is compared with the corresponding theoretical extended afterglow light 
curve we compute (dotted line). Also in this case we have a good correspondence between data and theoretical results. For 
completeness, the solid line shows again the theoretical extended afterglow light curve in the 15-150 keV energy range presented 
in Fig. 20. 



OPEN ISSUES IN CURRENT THEORETICAL MODELS 

The "fireshell" model addresses mainly the y and X-ray emission, which are energetically the most relevant part of the 
GRB phenomenon. The model allows a detailed identification of the fundamental three parameters of the GRB source: 
the total energy, the baryon loading, as well as the CBM density, filamentary structure and porosity. The knowledge 
of these phenomena characterizes the region surrounding the black hole up to a distance which in this source reaches 
~ 10 17 -10 18 cm. When applied, however, to larger distances, which corresponds to the latest phases of the X-ray 
afterglow, since the beginning of the "plateau" phase, the model reveals a different regime which has not yet been 
fully interpreted in its astrophysical implications. To fit the light curve in the soft X-ray regime for r>4x 10 17 cm, 
we must appeal to an enhancement of about six orders of magnitude in the g% factor (see above, and Fig. 21). This 
would correspond to a more diffuse CBM structure, with a smaller porosity, interacting with the fireshell. This points 
to a different leading physical process during the latest X-ray afterglow phases. When we turn to the optical, IR and 
radio emission, the fireshell model leads to a much smaller flux than the observed one, especially for r ~ 10 17 -10 18 
cm. Although the optical, IR, and radio luminosities have a minority energetic role, they may lead to the identification 
of crucial parameters and new phenomena occurring in the source, and they deserve maximum attention. 

In these latest phases for r > 10 17 cm it is currently applied the treatment based on synchrotron emission pioneered 
by Meszaros and Rees [124] even before the discovery of the afterglow [74]. Such a model has been further developed 
[see Refs. 125, 27, 28, and references therein]. Also in this case, however, some difficulties remain since it is necessary 
to invoke the presence of an unidentified energy injection mechanism [126]. Such a model appears to be quite 
successful in explaining the late phases of the X-ray emission of GRB 060614, as well as the corresponding optical 
emission, in terms of different power-law indexes for the different parts of the afterglow light curves [109, 127]. 
However, also in this case an unidentified energy injection mechanism between ~ 0.01 days and ~ 0.26 days appears 
to be necessary [127]. 

The attempt to describe the prompt emission via the synchrotron process by the internal shock scenario [see e.g. 
Refs. 52, 27, 28, and references therein] also encounters difficulties: Kumar and McMahon [53] have shown that the 
traditional synchrotron model can be applied to the prompt emission only if it occurs at r > 10 17 cm. A proposed 
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Figure 23. Energy emitted in: the P-GRB (solid line) and extended afterglow (dashed line), in function of the B parameter. 
The dotted line is the sum of the two lines. The two vertical lines represents the values of the B parameter of our two analyses, 
respectively. In our second analysis GRB 050509B results to be a "genuine" short. 



way-out of this problem, via the inverse Compton process, suffers of an "energy crisis" [see e.g. Ref. 128]. 

Interestingly, the declared region of validity of the traditional synchrotron model (r > 10 17 cm) is complementary 
to the one successfully described by our model (r < 10 17 -10 18 cm). Astrophysically, Xu et al. [127] have reached, 
within the framework of the traditional synchrotron model, two conclusions which are consistent with the results of 
our analysis of GRB 060614. First, they also infer from their numerical fit a very low density environment, namely 
n cbm ~ 0.04 particles/cm 3 . Second, they also mention the possibility that the progenitor of GRB 0606014 is a merging 
binary system formed by two compact objects. 

THE SEARCH FOR A "GENUINE SHORT" GRB: THE CASE OF GRB 050509B 

As we already discussed above, within the fireshell model the baryon loading is the key parameter to classify GRBs: 
if B < 10~ 5 we have what we call "genuine" short GRBs. In order to investigate if this is indeed the case of GRB 
050509B [129] we performed two different analyses, respectively with B = 3.7 x 10~ 3 and with B = 1.1 x 10~ 4 (see 
Fig. 23). 

Analysis 1. We identify the prompt emission of this GRB [see Ref. 129] with our P-GRB. Consequently, the 
extended afterglow corresponds to the observed X-ray afterglow (see Fig. 24). In this case, we have the total energy 
of the GRB estimated in Ef^ = 2.11 x 10 48 erg (which is a low energy for GRBs) and the baryon loading is 
B = 3.7 x 10~ 3 . With this choice of the fireshell parameters, we obtain that the P-GRB energy is Ep—grb — 1.6 x 10 4/ 
erg. More than 90% of the total energy is released in the extended afterglow, hence GRB 050509B cannot be classified 
as "genuine" short GRB. The reason for the non observability of the peak of this extended afterglow is that it results 
under the BAT threshold in the gamma-ray energy band, and before the beginning of the XRT observations (which 
sets around 100 s). 

Analysis 2. We performed an alternative analysis interpreting GRB 050509B as a "genuine" short GRB. We fit the 
BAT observations as the peak of the extended afterglow (see Fig. 25). In this case the total energy is E* ot — 5.07 x 10 49 
erg, the baryon loading is B = 1.1 x 10~ 4 , and this implies that the energy emitted on the P-GRB is almost 60% of the 
total one, Ep-grb = 3.30 x 10 49 erg. Differently with the previous case, it was not observable by BAT since its peak 
energy would be about 817 keV. According this second interpretation, GRB 050509B is a "genuine" short GRB. 
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Figure 24. Analysis 1 : the P-GRB corresponds to the BAT observations and the extended afterglow, that has a total energy that is 
much greater than the P-GRB one, to the XRT observations. 
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Figure 25. Analysis 2: both the BAT observations and the XRT ones are identified with the extended afterglow emission. The 
P-GRB flux is more then twice the extended afterglow one, but it is too hard to be observed by BAT. 



GRB 050509B within the Amati relation. In order to discriminate between the above analyses we checked if their 
results are compatible with the Amati relation [13, 20]. According to the fireshell model, only the extended afterglow 
component should satisfy the Amati relation, while the P-GRB component should not. This accounts for the fact that 
the "short" GRBs are outliers of the correlation. Therefore, according to the first interpretation the P-GRB which 
coincides with the BAT observation should be out from the correlation: this is indeed true. 

On the other hand, the second analysis reveals that the BAT observation should satisfy the correlation. Hence, this 
possibility have to be ruled out. We continue investigating the first possibility in order to obtain further information 
that from the astrophysical setting of this GRB constraint better the fit of the extended afterglow. 

Apart from this result, the new generation of high energy satellites are very important for the observation and study 
of the P-GRB, and for the identification of 'genuine" short GRBs. 




Figure 26. The Amati relation with our predictions about GRB 050509B in the two analyses. In the first one, the P-GRB results out 
from the correlation, as it should be. In the second one, the peak of the extended afterglow emission should satisfy the correlation, 
but it does not happen. Hence, the first analysis turns to be more correct. 



GRBS AND SNE: THE INDUCED GRAVITATIONAL COLLAPSE 

The Collapsar model [93, 130, 131, 132] proposes that GRBs arise from the collapse of a single Wolf-Rayet star 
endowed with fast rotation. This idea is purported by the evidence that many GRBs are close to star-forming regions 
and that this suggests that GRBs are linked to cataclysmic deaths of massive stars (M > 30M©). In such a model very 
massive stars are able to fuse material in their centers all the way to iron, at which point a star cannot continue to 
generate energy by fusion and collapses, in this case, immediately forming a black hole. Matter from the star around 
the core rains down toward the center and (for rapidly rotating stars) swirls into a high-density accretion disk. The 
mass of the accretion disk is around 0. 1 M©. The infall of this material into the black hole is assumed to drive a pair of 
jets (with opening angles < 10 degrees) out along the rotational axis, where the matter density is much lower than in 
the accretion disk, toward the poles of the star at velocities approaching the speed of light, creating a relativistic shock 
wave at the front [46]. The processes of core collapse and of accretion along the polar column and the jet propagation 
through the stellar envelope all together last <~ 10 sec [131]. The jet, as it passes through the star, is modulated by its 
interaction with the surrounding medium. In this way the Collapsar model attempts to explain the time structure of 
GRB prompt emission and to produce the variable Lorentz factor necessary for the internal shocks occurrence [112]. 
Moreover it is a prediction of this model that the central engine remains active for a long time after the principal burst 
is over, potentially contributing to the GRB afterglow [133]. 

Three very special conditions are required for a star to evolve all the way to a gamma-ray burst according to this 
theory: the star should be very massive (25M© Woosley [93], 35 — 40M© on the main sequence Fryer et al. [134]) to 
form a central black hole, the star rapidly rotates to develop an accretion torus capable of launching jets, and the star 
should have low metallicity in order to strip off its hydrogen envelope so the jets can reach the surface. As a result, 
gamma-ray bursts are far rarer than ordinary core-collapse supernovae, which only require the star to be massive 
enough to fuse all the way to iron. 

The consensus and the difficulties for the Collapsar model can be simply summarized: 

• Langer [135] asserts that long gamma-ray bursts are found in systems with abundant recent star formation, low 
metallicity environment. 

• The second evidence in favor of the Collapsar model is that there are several observed cases where a supernova 
is practical coeval with GRBs. 

• However, strong evidence against the Collapsar model comes from the fact that there were recently discovered two 
nearby long gamma-ray bursts which lacked a signature of any type Ib/c supernova: both GRB 060614 [88, 1 10, 
see also above] and GRB 060505 [1 14] defied predictions that a supernova would emerge despite intense scrutiny 




Figure 27. A sketch summarizing the induced gravitational collapse scenario. 



from ground-based telescopes. 

Within our fireshell model, we recall, the approach is drastically different, as already introduced in Ruffini et al. [6]. 
In fact in this framework, the SN which is often observed in temporal and spatial coincidence with the GRB cannot 
be interpreted as its progenitor because of the high quantity of ejected matter from the supernova explosion would 
prevent the GRB occurrence. Moreover: 

• It is very unlikely that a core collapse SN produces directly a black hole. 

• GRBs originate from gravitational collapses to black holes (see above). The possible explanation for the GRB-SN 
connection proposed in Ruffini et al. [6] is that both the GRB and the supernova progenitors belong to a binary 
system. Under special conditions it is possible that the GRB emission triggers the supernova explosion of the 
companion star. Alternatively, it is possible that the process of gravitational collapse to a black hole producing 
the GRB is "induced" by the supernova Ib/c on a companion neutron star [see Fig. 27 and Refs. 108, 51]. The 
faintness of this GRB class could be in this case naturally explained by the formation of the smallest possible 
black hole, just over the critical mass of the neutron star [136]. Moreover these systems occur in a low density 
CBM (1(T 2 -1 particle/cm 3 ). 

• Also the observation of the occurrence of "long" GRBs in star forming regions is explained by identifying the 
progenitor with a binary system formed by a neutron star and a star evolved out of the main sequence. 



Application to GRB 060218 



GRB 060218 triggered the BAT instrument of Swift on 18 February 2006 at 03:36:02 UT and has a T 9Q = 
(2100 ± 100) s [137]. The XRT instrument [138, 137] began observations ~ 153 s after the BAT trigger and continued 
for <~ 12.3 days [139]. The source is characterized by a fiat y-ray light curve and a soft spectrum [140]. It has an X-ray 
light curve with a long, slow rise and gradual decline and it is considered an X-Ray Flash (XRF) since its peak energy 
occurs at E p — 4.9+q 4 keV [141]. It has been observed by the Chandra satellite on February 26.78 and March 7.55 UT 
(t ~ 8.8 and 17.4 days) for 20 and 30 ks respectively [142]. The spectroscopic redshift has been found to be z = 0.033 
[143, 144]. The corresponding isotropic equivalent energy is Ei so = (1.9 ±0.1) x 10 49 erg [139] which sets this GRB 
as a low luminous one, consistent with most of the GRBs associated with SNe [106, 107, 145]. 

GRB 060218 is associated with SN2006aj whose expansion velocity is v ~ 0.1c [102, 146, 147, 107]. The host 
galaxy of SN2006aj is a low luminosity, metal poor star forming dwarf galaxy [148] with an irregular morphology 
[149], similar to the ones of other GRBs associated with SNe [150, 143]. 



The fit of the observed data 

In this section we present the fit of our fireshell model to the observed data (see Figs. 28, 31). The fit leads to a 
total energy of the e^ plasma E'J" = 2.32 x 10 50 erg, with an initial temperature T = 1 .86 MeV and a total number of 
pairs N e ± = 1 .79 x 10 55 . The second parameter of the theory, B = 1 .0 x 10~ 2 , is the highest value ever observed and is 
close to the limit for the stability of the adiabatic optically thick acceleration phase of the fireshell [for further details 
see Ref. 30]. The Lorentz gamma factor obtained solving the fireshell equations of motion [42, 45] is y = 99.2 at the 
beginning of the extended afterglow phase at a distance from the progenitor r = 7.82 x 10 12 cm. It is much larger 
than 7^5 estimated by Kaneko et al. [151] and Toma et al. [152]. 

In Fig. 28 we show the extended afterglow light curves fitting the prompt emission both in the BAT (15-150 keV) 
and in the XRT (0.3-10 keV) energy ranges, as expected in our "canonical GRB" scenario [51]. Initially the two 
luminosities are comparable to each other, but for a detector arrival time t d a > 1000 s the XRT curves becomes 
dominant. The displacement between the peaks of these two light curves leads to a theoretically estimated spectral 
lag greater than 500 s in perfect agreement with the observations [153]. We obtain that the bolometric luminosity in 
this early part coincides with the sum of the BAT and XRT light curves (see Fig. 31) and the luminosity in the other 
energy ranges is negligible. 

We recall that at t% <~ 10 4 s there is a sudden enhancement in the radio luminosity and there is an optical luminosity 
dominated by the SN2006aj emission [141, 142, 154]. Although our analysis addresses only the BAT and XRT 
observations, for r > 10 18 cm corresponding to t% > 10 4 s the fit of the XRT data implies two new features: 1) a 
sudden increase of the 8% factor from 3% = 1.0 x 10~ n to 8% = 1.6 x 10~ 6 , corresponding to a significantly more 
homogeneous effective CBM distribution (see Fig. 32b); 2) an XRT luminosity much smaller than the bolometric one 
(see Fig. 31). These theoretical predictions may account for the energetics of the enhancement of the radio and possibly 
optical and UV luminosities. Therefore, we identify two different regimes in the extended afterglow, one for^< 10 4 s 
and the other for t% > 10 4 s. Nevertheless, there is a unifying feature: the determined effective CBM density decreases 
with the distance r monotonically and continuously through both these two regimes from Yi c \y m — 1 particle/cm 3 at 
r = r a to n chm = 10~ 6 particle/cm 3 at r = 6.0 x 10 18 cm: n cbm « r~ a , with 1.0 < a < 1.7 (see Fig. 32a). 

Our assumption of spherical symmetry is supported by the observations which set for GRB 060218 an opening 
beaming angle larger than ~ 37° [106, 141, 142, 145]. 



The procedure of the fit 

The arrival time of each photon at the detector depends on the entire previous history of the fireshell [4]. Moreover, 
all the observables depends on the EQTS [47, 42] which, in turn, depend crucially on the equations of motion of the 
fireshell. The CBM engulfment has to be computed self-consistently through the entire dynamical evolution of the 
fireshell and not separately at each point. Any change in the CBM distribution strongly influences the entire dynamical 
evolution of the fireshell and, due to the EQTS structure, produces observable effects up to a much later time. For 
example if we change the density mask at a certain distance from the black hole we modify the shape of the lightcurve 
and consequently the evolution changes at larger radii corresponding to later times. Anyway the change of the density 
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Figure 28. GRB 060218 prompt emission: a) our theoretical fit (dotted line) of the BAT observations in the 15-150 keV energy 
band (dotted points); b) our theoretical fit (solid line) of the XRT observations in the 0.3-10 keV energy band (solid points) [Data 
from Ref. 141]. 



is not the only problem to face in the fitting of the source, in fact first of all we have to choose the energy in order 
to have Lorentz gamma factor sufficiently high to fit the entire GRB. In order to show the sensitivity of the fitting 
procedure I also present two examples of fits with the same value of B and different value of E'± . 

The first example has an E'J± = 1 .36 x 10 50 erg . This fit resulted unsuccessfully as we see from the Fig. 29, because 
the bolometric lightcurve is under the XRT peak of the extended afterglow. This means that the value of the energy 
chosen is too small to fit any data points after the peak of the extended afterglow. So we have to increase the value of 
the Energy to a have a better fit. In fact the parameters values have been found with various attempt in order to obtain 
the best fit. 

The second example is characterized by E'± = 1 .61 x 10 50 erg and the all the data are fitted except for the last point 
from 2.0 x 10 2 s to the end (see Fig. 30). I attempt to fit these last points trying to diminuishes the R values in order to 
enhance the energy emission, but again the low value of the Lorentz gamma factor, that in this case is 3 prevent the 
fireshell to expand. So again in this case the value of the Energy chosen is too small, but it is better than the previous 
attempt. In this case we increased the energy value of the 24%, but it is not enough so we decide to increase 16%. 
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Figure 29. GRB 060218 light curves with E'fJ = 1.36 x 10 50 erg: our theoretical fit (dotted line) of the 15-150 keV BAT 
observations and our theoretical fit (solid line) of the 0.3-10 keV XRT observations are represented (Data from: Campana et al. 
[141], Soderbergetal. [142]). 



So the final fit is characterized by the B = 1 .0 x 10~ 2 and by the E"£ = 2.32 x 10 50 erg. With this value of the energy 
we are able to fit all the experimental points. 



The fireshell fragmentation 

GRB 060218 presents different peculiarities: the extremely long Tw, the very low effective CBM density decreasing 
with the distance and the largest possible value of B = 10~ 2 . These peculiarities appear to be correlated. Following 
Ruffini et al. [108], we propose that in the present case the fireshell is fragmented. This implies that the surface of 
the fireshell does not increase any longer as r 2 but as r$ with j3 < 2. Consequently, the effective CBM density n c b m is 
linked to the actual one n a c f m by: 

n cb m = ^shelin a c c b t m , with ,$> she ii = (r*/r) a , (12) 

where r* is the starting radius at which the fragmentation occurs and a — 2 — j3 (see Fig. 32a). For r* — r Q we have 
Kbm = 1 particles/cm 3 , as expected for a "canonical GRB" [3] and in agreement with the apparent absence of a massive 
stellar wind in the CBM [142, 154, 155]. 

The & parameter defined in Eq.(13) has to take into account both the effect of the fireshell fragmentation {£% s heii) 
and of the effective CBM porosity (& c bm)'- 

3>,=M sheU xM cbm . (13) 
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Figure 30. GRB 060218 light curves with E'fJ = 1.61 x 10 50 erg: our theoretical fit (dotted line) of the 15-150 keV BAT 
observations and our theoretical fit (dashed line) of the 0.3-10 keV XRT observations are represented together with our theoretically 
computed bolometric luminosity (solid line). Data from: Campana et al. [141], Soderberg et al. [142]. 



The phenomenon of the dumpiness of the ejecta, whose measure is the filling factor, is an aspect well known in 
astrophysics. For example, in the case of Novae the filling factor has been measured to be in the range 10 2 — 10 5 
[156]. Such a filling factor coincides, in our case, with ^shell- 



Binaries as progenitors of GRB- SN systems 

The majority of the existing models in the literature appeal to a single astrophysical phenomenon to explain both the 
GRB and the SN ["collapsar", see e.g. Ref. 1 12]. On the contrary, a distinguishing feature of our theoretical approach 
is to differentiate between the SN and the GRB process (see above). The GRB is assumed to occur during the formation 
process of a black hole. The SN is assumed to lead to the formation of a neutron star (NS) or to a complete disruptive 
explosion without remnants and, in no way, to the formation of a black hole. In the case of SN2006aj the formation of 
such a NS has been actually inferred by Maeda et al. [104] because of the large amount of 58 Ni (0.05M©). Moreover 
the significantly small initial mass of the SN progenitor star M w 2OM is expected to form a NS rather than a black 
hole when its core collapses [104, 148, 157, 158]. In order to fulfill both the above requirement, we assume that the 
progenitor of the GRB and the SN consists of a binary system formed by a NS close to its critical mass collapsing to 
a black hole, and a companion star evolved out of the main sequence originating the SN. The temporal coincidence 
between the GRB and the SN phenomenon is explained in term of the concept of "induced" gravitational collapse 
[6, 108]. There is also the distinct possibility of observing the young born NS out of the SN [see e.g. Ref. 108, and 
references therein]. 

It has been often proposed that GRBs associated with SNe Ib/c, at smaller redshift 0.0085 < z < 0.168 [see e.g. 




Figure 31. GRB 060218 complete light curves: our theoretical fit (dotted line) of the 15-150 keV BAT observations, our 
theoretical fit (solid line) of the 0.3-10 keV XRT observations and the 0.3-10 keV Chandra observations are represented together 
with our theoretically computed bolometric luminosity (double dotted line) [Data from Refs. 141, 142], 



Ref. 105, and references therein], form a different class, less luminous and possibly much more numerous than the 
high luminosity GRBs at higher redshift [102, 103, 104, 105]. Therefore they have been proposed to originate from a 
separate class of progenitors [106, 107]. In our model this is explained by the nature of the progenitor system leading 
to the formation of the black hole with the smallest possible mass: the one formed by the collapse of a just overcritical 
NS [159, 108]. 



Conclusions 

GRB 060218 presents a variety of peculiarities, including its extremely large Tgo and its classification as an XRF. 
Nevertheless, a crucial point of our analysis is that we have successfully applied to this source our "canonical GRB" 
scenario. 

Within our model there is no need for inserting GRB 060218 in a new class of GRBs, such as the XRFs, alternative 
to the "canonical" ones. This same point recently received strong observational support in the case of GRB 060218 
[153] and a consensus by other models in the literature [151]. 

The anomalously long T^q led us to infer a monotonic decrease in the CBM effective density giving the first clear 
evidence for a fragmentation in the fireshell. This phenomenon appears to be essential in understanding the features of 
also other GRBs [see e.g. GRB 050315 in Refs. 108, 33]. 

Our "canonical GRB" scenario originates from the gravitational collapse to a black hole and is now confirmed over 
a 10 6 range in energy [see e.g. Ref. 3, and references therein]. It is clear that, although the process of gravitational 
collapse is unique, there is a large variety of progenitors which may lead to the formation of black holes, each one with 
precise signatures in the energetics. The low energetics of the class of GRBs associated with SNe, and the necessity of 
the occurrence of the SN, naturally leads in our model to identify their progenitors with the formation of the smallest 
possible black hole originating from a NS overcoming his critical mass in a binary system. For GRB 060218 there is 
no need within our model for a new or unidentified source such as a magnetar or a collapsar. 




Distance (cm) 



Figure 32. The CBM distribution parameters: a) the effective CBM number density (solid line) monotonically decreases with the 
distance r following Eq.(12) (dashed line); b) the 3% parameter vs. distance. 



GRB 060218 is the first GRB associated with SN with complete coverage of data from the onset all the way up to 
~ 10 6 s. This fact offers an unprecedented opportunity to verify theoretical models on such a GRB class. For example, 
GRB 060218 fulfills the Amati et al. [13] relation unlike other sources in its same class. This is particularly significant, 
since GRB 060218 is the only source in such a class to have an excellent data coverage without gaps. We are currently 
examining if the missing data in the other sources of such a class may have a prominent role in their non-fulfillment 
of the Amati et al. [13] relation [Dainotti et al., in preparation; see also Ref. 160]. 



THEORETICAL BACKGROUND FOR GRBS' EMPIRICAL CORRELATIONS 

The detection of GRBs up to very high redshifts [up to z = 6.7, see Ref. 10], their high observed rate of one every few 
days, and the progress in the theoretical understanding of these sources all make them useful as cosmological tools, 
complementary to supernovae la, which are observed only up to z = 1.7 [11, 12]. One of the hottest topics on GRBs 
is the possible existence of empirical relations between GRB observables [13, 14, 15, 16, 17, 18], which may lead, if 



confirmed, to using GRBs as tracers of models of universe. The first empirical relation, discovered when analyzing the 
BeppoSAX so-called "long" bursts with known redshift, was the "Amati relation" [13]. It was found that the isotropic - 
equivalent radiated energy of the prompt emission Ei so is correlated with the cosmological rest-frame vFy spectrum 
peak energy E p /. E p j °= with a = 0.52 ± 0.06 [13]. The existence of the Amati relation has been confirmed by 

studying a sample of GRBs discovered by Swift, with a = 0.49+Q Q5 [19, 20]. 

Swift has for the first time made it possible to obtain high quality data in selected energy bands from the GRB 
trigger time all the way to the latest extended afterglow phases [161]. This has given us the opportunity to apply our 
theoretical "fireshell" model, thereby obtaining detailed values for its two free parameters, namely for the total energy 
Ef ot and the baryon loading B of the fireshell, as well as for the effective density and filamentary structure of the 
CBM. From this we were able to compute multi-band light curves and spectra, both instantaneous and time-integrated, 
compared with selected GRB sources, such as GRB 050315. 

In the "fireshell" model, Ef ot comprises two different components: (i) the P-GRB with energy E P _ GRB , emitted at 
the moment when the e + e~ -driven accelerating baryonic matter reaches transparency, and (ii) the following extended 
afterglow phase with energy E a f t , with the decelerating baryons interacting with the CBM [5]. These two phases are 
clearly distinguishable by their relative intensity and temporal separation in arrival time. We have 

E tot =Ep- GRB + Eaft ■ (14) 

What is usually called the "prompt emission" corresponds within the fireshell model to the P-GRB together with the 
peak of the extended afterglow [see below, e.g. Ref. 5, 71, 3, 51, 33, 162, 34, and references therein]. 

Among the crucial issues raised by the Amati relation, there are its theoretical explanation and its possible depen- 
dence on the assumed cosmological model. We examined a set of "gedanken" GRBs, all at the same cosmological 
redshift of GRB 050315. Such a set assumes the same fireshell baryon loading and effective CBM distribution as 
GRB 050315 [71] and each "gedanken" GRB differs from the others uniquely by the value of its total energy E e m . We 
then considered a second set of "gedanken" GRBs, differing from the previous one by assuming a constant effective 
CBM density instead of the one inferred for GRB 050315. In both these sets, we looked for a relation between the 
isotropic-equivalent radiated energy of the entire extended afterglow E a f t and the corresponding time-integrated vF v 
spectrum peak energy E p : 

E p oc{E af1 f. (15) 

In this chapter, after briefly recalling the various spectral-energy correlations mentioned above, we present the 
derivation of the E p - E a f t relation for the two sets of "gedanken" GRBs. 



Spectral-energy correlations 

Many empirical spectral-energy correlations exist, some are purely phenomenological and assumption free while 
others are based on assumptions and are dependent on model, basically the standard fireball model [43]. Some 
correlations assume spherical symmetry while others assume collimated (jet) emission. This last case was triggered 
by the observation by Frail et al. [163] that the collimation corrected energetics of those GRBs of know jet aperture 
angles clustered into a narrow distribution, Ey = (1 — cos 6j)Ei s0 <~ 10 51 erg. The opening angle of the jet is estimated 
within the standard model as 

e } = o.ieil^Y'V^] ; h 

a = 0.2016 1^-) (-11—) ; W (16) 

where fjeu is the break time measured in days and z is the redshift. The efficiency T]y relates the isotropic kinetic 
energy of the fireball £k,iso to the prompt emitted energy E{ s0 : £k,iso = E^ /r]y. Usually, it is assumed a constant value 
for all bursts, i.e. rjy = 0.2 (after its first use by Frail et al. [163], following the estimate of this parameter in GRB 
970508 [164]). In the homogeneous (H) case, n is the CircumBurst density, independent from the radial coordinate; 
for the wind (W) case, the density is a function of the radial coordinate, n(r) =Ar~ 2 and A* is the value of A 
(A = M w /(4nv w ) = 5 x 10 n AH, g cm" 1 ) when setting the wind mass loss rate to M w = 10~ 5 M Q ;yr~ 1 and the wind 
velocity to v w = 10 3 km s _1 . Usually, a constant value (i.e. A* = 1) is adopted for all bursts. 



l+zj \E is0 ,52 

t- a \ l ' A / 71 A \ l l A 
'jet.d \ / HyA* 



The most important spectral-energy correlations are: 

• The Amati relation: It was historically the first correlation discovered, considering BeppoSAX bursts [ 1 3] . It was 
found that the isotropic-equivalent radiated energy of the prompt emission Ei so is correlated with the cosmological 
rest-frame vF v spectrum peak energy E p /. E p j « EfJ. This correlation, recently updated [20] to a larger sample, 
holds for all but two long bursts, while no short burst satisfies it. The long burst outliers are GRB 980425 and 
the debated GRB 031203 [165]. As far as short bursts are concerned, there are two cases: the burst with an initial 
spike-like emission followed by a soft bump (short burst with afterglow) and the short burst with no afterglow. 
The burst belonging to the first class are what we named [33] "fake" or "disguised" short GRBs, while the ones 
belonging to the second case are the "genuine" short GRBs. Both classes, as already said above, does not follow 
the Amati relation, but, if one excludes the initial spike-like emission and considers only the soft later part of the 
bursts in the first class, then the Amati relation is recovered [83, 81]. 

• The Yonetoku correlation: Yonetoku et al. [15] showed that also the peak luminosity L p j S0 of the prompt 

emission correlates with E p , in the same way as Ei so '. Ep oc E • . The scatter is similar to the scatter of the 
Amati correlation, and the outliers are the same as well. 

• The Ghirlanda correlation: Assuming a collimated emission, Ghirlanda et al. [14] found that the collimation 
corrected (by a factor (1 — cos0 7 )) energy, E y , is tightly correlated with E p . The correlation is E p j °< Ey 7 . As 
outlined above, this relation is based on a theoretical model needed to calculate dj, that in turns relies on the 
assumptions adopted for the efficiency and the CircumBurst density and profile. 

• The Liang & Zhang correlation: To find the jet angle 9j, as explained above, one needs a model and some 
assumptions; the Liang and Zhang [16] correlation instead is entirely phenomenological, so model independent 
and assumptions free. It involves three observables (plus the redshift) and it is of the form E iso °c Ep jj . It is 
consistent [166] with the Ghirlanda correlation, and has similar spread. 

• The Firmani correlation: The Firmani et al. [17] correlation links three quantities of the prompt emission: 
the bolometric isotropic peak luminosity L p , the peak energy E p i so of the time integrated spectrum, and a 
characteristic time: 7o.45, which is the time interval spanned by the brightest 45% of the total light curve counts 
above the background. This time is used to characterize the variability properties of the prompt emission [167]. 

3/2 —1/2 

The correlation is of the form: L p ,i so E p ' T Q 45 ' . Also this relation is model independent and assumption free. 

The E p - E a f t relation 

In our approach, only the entire extended afterglow emission is considered in establishing our E p - E a f, relation. 
From this assumption one derives, in a natural way, that the Amati relation holds only for long GRBs, where the 
P-GRB is negligible, and not for short GRBs [115]. 

We can compute the "instantaneous" spectrum of GRB 050315 at each value of the detector arrival time during the 
entire extended afterglow emission. Such a spectrum sharply evolves in the arrival time, presenting a typical hard- 
to-soft behavior [71]. We then computed the vF v time-integrated spectrum over the total duration of our extended 
afterglow phase, that is, from the end of the P-GRB up to when the fireshell reaches a Lorentz gamma factor close to 
unity. We can then define the energy E p as the energy of the peak of this vF v time-integrated spectrum, and we look 
at its relation with the total energy E a f, of the extended afterglow. 

We construct two sets of "gedanken" GRBs at a fixed cosmological redshift, therefore independently of the 
cosmological model. The first set assumes the same fireshell baryon loading and effective CBM distribution as GRB 
050315 (see Fig. 33) and each "gedanken" GRB differs from the others uniquely by the value of its total energy 
Ef ot . The second set assumes a constant effective CBM density <~ 1 particle/cnr instead of the one inferred for GRB 
050315. 

In our model, E a f, is a fixed value determined by Ef ot and B, so clearly there are no errors associated to it. 
Instead, E p is evaluated from the numerically calculated spectrum, and its determination is therefore affected by 
the numerical resolution. Choosing a 5% error on E p , which is consistent with our numerical resolution, we checked 
that this value is reasonable looking at each spectrum. Figure 34 shows the time-integrated spectrum corresponding to 
Ef ot = 3.40 x 10 51 erg with the error around E p . 
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Figure 33. The effective CBM number density inferred from the theoretical analysis of GRB 050315. Details in Ruffini et al. 
[71]. 
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Figure 34. The vF v time-integrated spectrum over the total duration of our extended afterglow phase for the "gedanken" GRB of 
the first set with total energy Ef ot = 3.40 x 10 51 . The two vertical lines constrain the 5% error region around the peak. We determine 
£ p = 5.82keV±5%. 

Results and discussion 

Figure 35 shows the E p - E a f, relation of the "gedanken" GRBs belonging to the first set (red points). It extends 
over two orders of magnitude in energy, from 10 51 to 10 53 erg, and is well-fitted by a power law E p °< (E a f t ) a 
with a = 0.45 ±0.01. We emphasize that such a power-law slope strictly agrees with the Amati relation, namely 
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Figure 35. The E p - E a j, relation: the results of the simulations of the first set of "gedanken" GRBs (points marked as crosses) 
are well-fitted by a power law (solid line) E p oc (E a f t ) a with a = 0.45 ± 0.01. The points marked as "X" are the results of the 
extension of the first set above 10 53 erg. 
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Figure 36. The vF v time-integrated spectrum over the total duration of our extended afterglow phase for the "gedanken" GRB of 
the extended first set with total energy Ef ot = 6.95 x 10 53 erg. The vertical lines constrain the 5% error region around each peak. 

E p i ex (E iso ) a , with a = 0.49+qqj [20]. We recall that E p is the observed peak energy; i.e., it is not rescaled for the 
cosmological redshift, because all the "gedanken" GRBs of the set are at the same redshift of GRB 050315, namely 
z = 1 .949 [100]. The normalization is clearly different from the Amati one. 

If we try to extend the first sample of "gedanken" GRBs below 10 51 erg, the relevant CBM distribution would be 
for r < 10 16 cm, where no data are available from the GRB 050315 observations. If we try to extend the first set 
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Figure 37. The second set of "gedanken" GRBs. Clearly, in this case there in no relation between E p and E a f t . 

of "gedanken" GRBs above 10 53 erg, we notice that for Ef ot > 10 54 erg the small "bump", which can be noticed 
between 0.2 and 1 .0 keV in the spectrum of Fig. 34, evolves into a low-energy second spectral peak that is even higher 
than the high-energy one (see Fig. 36). We are currently investigating whether this low-energy second peak is a real, 
theoretically predicted spectral feature that may be observed in the future in highly energetic sources. There is also 
the other possibility that the low-energy and late part of our GRB 050315 fit is not enough constrained by the XRT 
observational data so that this effect is magnified by the Ef at rescaling. 

The high-energy spectral peak is due to the emission at the peak of the extended afterglow, and therefore due to the 
so-called "prompt emission". The low-energy one is due to late-time soft X-ray emission. Therefore, the high-energy 
spectral peak is the relevant one for the Amati relation. We find indeed that such a high-energy spectral peak still 
fulfills the E p - E a r t relation for Ef ot ~ 10 54 erg, with a possible saturation for Ef ot > 10 54 erg (see Fig. 35). 

Figure 37 clearly shows that in the second set of "gedanken" GRBs, built assuming a constant effective CBM density 
~ 1 particle/cm 3 , instead of the one specifically inferred for GRB 050315, there in no relation between E p and E a f t . 

Conclusions 

The high-quality Swift data, for the first time giving gapless and multiwavelength coverage from the GRB trigger all 
the way to the latest extended afterglow phases, have led to a complete fit of the GRB 050315 multiband light curves 
based on our fireshell model. We fixed the free parameters describing the source and determined the instantaneous and 
time-integrated spectra during the entire extended afterglow. 

Starting from this, we examined two sets of "gedanken" GRBs, constructed at a fixed cosmological redshift. The first 
set assumes the same fireshell baryon loading and effective CBM distribution as GRB 050315, and each "gedanken" 
GRB differs from the others uniquely by the value of its total energy Ef ot . The second set assumes a constant effective 
CBM density ~ 1 particle/cm 3 instead of the one inferred for GRB 050315. 

Recalling that the "canonical" GRB light curve in the fireshell model is composed of two well-separated compo- 
nents, the P-GRB, and the entire extended afterglow, we looked for a relation in both sets between the isotropic- 
equivalent radiated energy of the entire extended afterglow E a f, and the corresponding time-integrated vF v spectrum 
peak energy E p : E p « (E a f t ) a . In doing so, we assumed that the Amati relation is directly linked to the interaction 
between the accelerated baryons and the CBM. The P-GRBs, which originate from the fireshell transparency, do not 
fulfill the Amati relation in our approach. Consequently, the short GRBs, which have a vanishing extended afterglow 



with respect to the P-GRB, should also not fulfill the Amati relation. This last point is supported by the observational 
evidence [20]. 

We notice that the first set of "gedanken" GRBs fulfills the E p °c {E a f t ) a relation very well with a = 0.45 ± 0.01. 
This slope strongly agrees with the Amati relation. In contrast, no relation between E p and E a f t seems to hold for the 
second set. We conclude that the Amati relation originates from the detailed structure of the effective CBM. 

Turning now to the analogies and differences between our E p - E a f t relation and the Amati one, our analysis 
excludes the P-GRB from the prompt emission, extends all the way to the latest extended afterglow phases, and is 
independent of the assumed cosmological model, since all "gedanken" GRBs are at the same redshift. The Amati 
relation, on the other hand, includes the P-GRB, focuses only on the prompt emission, being therefore influenced by 
the instrumental threshold that fixes the end of the prompt emission, and depends on the assumed cosmology. This 
might explain the intrinsic scatter observed in the Amati relation [20]. Our theoretical work is a first unavoidable step 
toward supporting the use of the empirical Amati relation for measuring the cosmological parameters. 

THERMALIZATION PROCESS OF ELECTRON-POSITRON PLASMA WITH 

BARYON LOADING 

Initial evolution of electron-positron-photon plasma in the source of a GRB has a key role in the subsequent dynamics 
of the fireshell. In particular, particle spectra, temperatures, chemical potentials all need to be known in order to 
describe acceleration of the fireshell and in general its expansion in terms of hydrodynamics. Since quite different 
theoretical arguments existed on the initial state of optically thick electron-positron-photon plasma in GRBs [see e.g. 
Refs. 168, 96] we turned to analysis of kinetic properties of nonequilibrium electron-positron pairs. 

Having this goal in mind Aksenov et al. [122] solved numerically relativistic Boltzmann equations for distribution 
functions of electrons, positrons and photons, assuming their uniform spatial distribution. Considering energy density 
in the range, typical for GRBs, the relevant thermalization timescales were determined. It turns out that particles reach 
kinetic equilibrium on a timescale t < 10~ 14 sec, when distribution functions of electrons/positrons (photons) acquire 
Fermi-Dirac (Bose-Einstein) form, all particles have a common temperature but nonzero chemical potentials. Further, 
on a timescale t < 10~ 12 sec chemical potentials vanish and particles reach thermal distribution. 

Since in many bursts baryon loading is dynamically significant Aksenov et al. [40] considered proton admixture 
parametrized by the parameter B = n p m p c 2 / p r , where n p is the number density of protons, m p is their mass, p r is the 
radiative energy density (including the energy density of electron-positron pairs). Independent on the baryon loading 
the thermalization timescale was found to be t < 10~ u sec for such a plasma, despite thermalization process is more 
complicated. 

As example, we show energy densities (Fig. 38), number densities (Fig. 39), temperatures and chemical potentials 
(Fig. 40, 41 respectively) depending on time. Initial and final spectra of particles are shown in Fig. 42. Initial conditions 
were chosen with flat spectral densities and total energy density p = 10 24 erg/cm 3 . This initial state is clearly far from 
equilibrium. Interactions between particles change distribution functions such that at the moment ti = 4 x 10~ 14 sec, 
shown by the vertical line on the left in Fig. 40 and 41, distribution functions acquire an equilibrium form with 
temperature of photons and pairs is 9^ — ksT^/ (m e c 2 ) ~ 1.5, while the chemical potentials of these particles are 
Vk = <Pk/(m e c 2 ) ~ —7, where k B is Boltzmann's constant, m e is electron mass, and c is the speed of light. These 
changes are essentially due to binary interactions, which do not change number of particles. Further evolution of 
the system is due to triple interactions. Only triple interactions, which do not conserve the number of particles, are 
able to change chemical potentials of particles. The chemical potential of photons reaches zero at = 10~ 12 sec, 
shown by the vertical line on the right in Fig. 40 and 41, which means that electrons, positrons and photons come to 
thermal equilibrium. Protons join thermal equilibrium the last. Simulations with different initial conditions show that 
thermalization timescale depends essentially on total energy density and the baryonic loading parameter B [for details 
seeRef. 169]. 

Such short timescales, compared to a typical expansion timescale t ex ~ Rq/c ~ 10~ 3 sec where Rq is initial size of 
the plasma allow to speak of completely thermalized plasma long before expansion starts [for details see Ref. 169]. 




Figure 38. Depencence on time of energy densities of electrons (green), positrons (red), photons (black) and protons (blue). Total 
energy density is shown by dotted black line. Interaction between pairs and photons operates on very short timescales up to 10 -23 
sec. Quasi-equilibrium state is established at ~ 10 -14 sec which corresponds to kinetic equilibrium for pairs and photons. Protons 
start to interact with them as late as at f t h ~ 10~ 13 sec. 



CRITICAL ELECTRIC FIELDS ON THE SURFACE OF MASSIVE CORES AND 
DYADOTORUS OF THE KERR-NEWMAN GEOMETRY 

Critical electric Fields on the surface of massive cores 

One of the most active field of research has been to analyse a general approach to Neutron Stars based on the 
Thomas-Fermi ultrarelativistic equations amply adopted in the study of superheavy nuclei. The aim is to have a 
unified approach both to superheavy nuclei, up to atomic numbers of the order of 10 s — 10 6 , and to what we have 
called "Massive Nuclear Cores", which are 

• characterized by atomic number of the order of 10 57 ; 

• composed by neutrons, protons and electrons in j3 -equilibrium; 

• expected to be kept at nuclear density by self gravity. 

The analysis of superheavy nuclei has historically represented a major field of research [170, 171, 172, 173, 174], 
guided by Prof. V. Popov and Prof. W. Greiner and their schools. This same problem was studied in the context of the 
relativistic Thomas-Fermi equation also by R. Ruffini and L. Stella [175, 176], already in the 80s. The recent numerical 
approach has shown the possibility to extrapolate this treatment of superheavy nuclei to the case of Massive Nuclear 
Cores [177]. The very unexpected result has been that also around these massive cores there is the distinct possibility 
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Figure 39. Depencence on time of concentrations of electrons (green), positrons(red), photons (black) and protons (blue). Total 
number density is shown by dotted black line. In this case kinetic equilibrium between electrons, positrons and photons is reached 
at ~ 1CP 14 sec. Protons join thermal equilibrium with other particles at f t j, ~ 4 x 1CP 12 sec. 



of having an electromagnetic field close to the critical value E c = ^fg-, although localized in a very narrow shell of 
the order of the electron Compton wavelength (see Fig. 43, 44). 

The welcome result has been that all the analytic work [178] developed by Prof. Popov and his Russian collaborators 
can be straightforwardly applied to the case of massive cores, and the j3 -equilibrium condition is properly taken into 
account. In Ref. [179], we show that globally neutral massive cores can be gravitationally bound and the value of 
the charge-to-mass ratios predicted at the surface of massive cores coincides with the range of values expected in 
astrophysical scenarios for Kerr-Newman black holes, in addition to a further verification of the over critical electric 
field on the massive core surface, numerically obtained by Ruffini, Xue and Rotondo already in 2007 [177]. A large 
variety of problems has emerged, and working progress in solving these problems have been going on in the direct 
discussion with or participation by Prof. Greiner, Prof. Popov, and Prof, 't Hooft at ICRANet center located at Pescara. 
The crucial issue to be debated is the stability of such cores under the competing effects of self gravity and Coulomb 
repulsion. In order to probe this stability, we have started a new approach to the problem within the framework of 
general relativity. The object of the work by Patricelli and Rueda is the generalization of the Tolman-Oppenheimer- 
Volkoff equation duly taking into account the elecrodynamical contribution. The major scientific issue here is to have a 
unified approach solving the coupled system of the general relativistic self gravitating electrodynamical problem with 
the corresponding formulation of the Thomas-Fermi equation in the framework of general relativity. Prof, 't Hooft, in 
a series of lectures at Pescara, has forcefully expressed the opinion that necessarily, during the process of gravitational 
collapse, it should occur a more extended distribution of the electromagnetic field to the entire core of the star and not 
only confined to a thin shell. This is a necessary condition in order to transmit the gravitational energy of the collapse 
to the electrodynamical component of the field giving possibly rise to large pair creation processes. This crucial idea 




Figure 40. Depencence on time of dimensionless temperature of electrons (green), positrons (red), photons (black) and protons 
(blue). The temperature for pairs and photons acquires physical meaning only in kinetic equilibrium at ~ 10~ 14 sec. Protons are 
cooled by the pair-photon plasma and acquire common temperature with it as late as at ~ 4 x 10 -12 sec. 



is currently being pursued by the application to this system of a classical work of Feynmann-Metropolis and Teller, 
who considered in relativistic Thomas-Fermi the crucial role of non-degeneracy. 



On the Dyadotorus of the Kerr- Newman Geometry 

In the merging process of two neutron stars and in the final process of gravitational collapse of a black hole it is 
possible that very large electromagnetic field strength larger than the critical value of vacuum polarization E c do occur 
[3]. The description of the time evolution of the gravitational collapse and the associated electrodynamical process 
(occurring on characteristic time scales T = GM/c 3 ~ 5 x 1O~ 5 M/M s) are too complex for a direct description. A 
more confined problem is the case of an already formed Kerr-Newman black hole. 

This deserves analysis in itself as a theoretical problem and may represent a physical condition asymptotically 
reached in the process of gravitational collapse. Such an asymptotic configuration will be reached when all the multi- 
poles departing from the Kerr-Newman geometry have been radiated away either by process of vacuum polarization 
or electromagnetic and gravitational waves. This simplified problem may lead to a direct evaluation of the energetics 
as well as of the created e~e + pairs occurring on time scales At = h/(m e c 2 ) ~ 10~ 21 s. 

Therefore we explore the initial condition for such a process by the definition of the spatial extent of a "dyadotorus" 
which generalizes to the Kerr-Newman geometry the concept of the "dyadosphere" previously introduced in the case 
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Figure 41. Depencence on time of dimensionless chemical potential of electrons (green), positrons (red), photons (black) and 
protons (blue). The chemical potential for pairs and photons acquires physical meaning only in kinetic equilibrium at ~ 10~ 14 
sec, while for protons this happens at fa — 4 x 10~ 12 sec. At this time chemical potential of photons has evolved to zero and thermal 
equilibrium has been already reached. 



of the spherically symmetric Reissner-Nordstrom geometry [180, 26]. 

Damour and Ruffini [23] showed that vacuum polarization processes a la Sauter-Heisenberg-Euler-Schwinger [181] 
can occur in the field of a Kerr-Newman black hole endowed with a mass ranging from the maximum critical mass for 
neutron stars (3.2M©) all the way up to 7.2 x 10 6 M Q . It is an almost perfectly reversible process in the sense defined 
by Christodoulou and Ruffini [182], leading to a very efficient mechanism of extracting energy from the black hole. 

In the case of absence of rotation in spacetime, we have a Reissner-Nordstrom black hole as the background 
geometry. The region where vacuum polarization processes take place is a sphere centered about the hole, and has 
been called dyadosphere [180, 26]. We investigate how the presence of rotation in spacetime modifies the shape of the 
surface containing the region where electron-positron pairs are created. 

Due to the axial symmetry we call that region as dyadotorus and we give the conditions for its existence. We have 
defined the dyadotorus as the locus of points where E = kE c with k some positive constant which can be less than one 
[see Ref. 183, for details]. We have found that the geometry of the dyadotorus is indeed torus-like when 
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k>- - T « 6.6 x 10 4 -^ , (17) 



where n = M/MQ, t, = Q/M and a = J/M 2 being M, Q and J the mass, the charge and angular momentum of the 
black hole. Otherwise, it becomes ellipsoid-like. This can be seen from Fig. 45. 

An estimate of the electromagnetic energy contained in the dyadotorus can be calculated by using, for example, three 
different definitions of it commonly adopted in the literature, i.e. the standard definition in terms of the timelike Killing 
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Figure 42. Spectral density as function of particle energy for electrons (green), positrons (red), photons (black) and protons (blue) 
in initial and final time moments of the computation. Fits of the spectra with chemical potentials and temperatures corresponding 
to thermal equilibrium state are also shown by yellow (electrons and positrons), gray (photons) and light blue (protons) thick lines. 
The final photon spectrum is black body one. 



vector [see e.g. Ref. 184], the one recently suggested by Katz, Lynden-Bell and Bicak [185, 186] for axially symmetric 
asymptotically flat spacetimes, which is an observer dependent definition of energy, and the last one involving the 
theory of pseudotensors [see e.g. 187]. All these approaches are shown to give the same results. 

From this, we find that in addition to the topological differences between the dyadotorus and the dyadosphere, larger 
field strengths and electromagnetic energy are allowed in the case of a Kerr-Newman geometry close to the horizon, 
when compared with a Reissner-Nordstrom black hole of the same mass energy and charge to mass ratio. 

SELECTED PROCESSES ORIGINATING HIGH-ENERGY EMISSION 

The knowledge of the radiation mechanisms is crucial for the correct understanding of many astronomical observa- 
tions. In particular in gamma ray astronomy the observational data can be often explained by two or more production 
mechanisms. It is therefore important to model correctly the different interactions producing gamma rays. Among all 
mechanisms producing gamma rays hadronic interactions between nucleons which produce pions, which in turn decay 
into photons and neutrinos, are one of the most studied models. High energy collisions of nucleons cannot be treated 
perturbatively because of the large value of the interaction constant in nuclear interactions. Many authors have occu- 
pied themselves with this problem and already in 1950 Fermi developed an elegant statistical method for computing 
the multiple production of particles in collisions of high energetic protons. In the meantime a very large set of data 
has been acquired from high energy accelerators. Our aim here is 1. to rederive the Fermi theoretical equations, 2. to 
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Figure 43. Number density of electrons, protons and neutrons. 
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Figure 44. Electric Field in units of the critical field. 





Figure 45. The projections of the dyadotorus on the X — Z plane corresponding to different values of the ratio |E|/£ C = k are 
shown in the left pane for /i = 10 and E, = 1.49 x 10~ 4 . The corresponding plot for the dyadosphere with the same mass energy M 
and charge to mass ratio £, is shown in the right pane for comparison. 



compare them to the experimental data and 3. to explore possibilities of observing such phenomena. 



Fermi's approach to the study of hadronic interactions 

In treating high energy collisions of nucleons, Fermi made the assumption that the possible final configurations 
of the system are determined by the statistical weights of the various possible final configurations and accordingly 
developed a statistical method to determine the final particles produced [21]. First of all one might think of many 
different final configurations for the system after the collision, but conservation laws of charge and of momentum, as 
well as the feasibly of the processes, have to be taken into account. Transitions in Yukawa's theory, in which charged 
and neutral pions are created, are therefore the most probable processes taking place. So during the collisions of high 
energetic hadrons a large amount of energy is released in a small volume around the hadrons and used to form pions. 
In view of the strong interactions between these pions, one can imagine that the energy available in the small volume 
will be rapidly distributed to the different pions having different energies. In other words the energy will be statistically 
distributed among all degrees of freedom of the system. Fermi himself said: "When two nucleons collide with very 
great energy in their center of mass system this energy will be suddenly released in a small volume surrounding the two 
nucleons. The event is a collision in which the nucleons with their surrounding retinue of pions hit against each other 
so that all the portion of space occupied by the nucleons and by their surrounding pion field will be suddenly loaded 
with a very great amount of energy. Since the interactions of the pion field are strong we may expect that rapidly this 
energy will be distributed among the various degrees of freedom present in this volume according to statistical laws. 
One can then compute statistically the probability that in this tiny volume a certain number of pions will be created 
with a given energy distribution. It is then assumed that the concentration of energy will rapidly dissolve and that the 
particles into which the energy has been converted will fly out in all directions." 

Fermi's method resembles somehow Heisenberg's approach [188] to treat high energy collisions of nucleons, with 
the difference that Heisenberg used qualitative ideas of turbulence whereas, Fermi believed that in high energetic 
processes statistical equilibrium is reached. 

According to Fermi the process proceeds as follows: in the laboratory frame a very energetic proton scatters off a 
proton target. For convenience the process is examined in the center of mass of the system. The only parameter which 
has to be tuned in Fermi's theory is the volume in which the energy is dumped. The value of this parameter can be 
modified to improve the agreement with the experiments. Fermi defined £2 being the volume at rest in the laboratory 
frame that contains the energy of the colliding particles. Since the particles mediating the Yukawa interactions are the 



pions, the volume is taken as a sphere with radius of order of the pion Compton wavelength A£ = h/m K c = 1 .4 x 10~ 13 
cm. 

We will use the subscript "0" to indicate quantities in the cm. frame, and no subscript to indicate the same 
quantities in the laboratory frame. Analogously, we use the superscript ' to indicate quantities after the collision, 
and no superscript to indicate the same quantities before the collision. 

Let's focus for the moment on the target proton, which is at rest in the laboratory frame. Its associated volume £2 is: 

&=\nR\ (18) 

with 

R = aX c K , (19) 

where a is a free parameter of the order of unity which Fermi leaves free to better fit the experimental data. The same 
volume measured in the cm. frame, where the target proton is moving with Lorentz factor 7, is given by: 

Q. = -CI. (20) 

r 

In order to determine the Lorentz factor 7, we consider the case in which, in the initial reference frame, the laboratory 
frame, a particle with mass m\ and energy E\ collides with a particle with mass m 2 which is at rest. The total energy 
of the two particles is: 

E =E l +E 2 =E l +m 2 c 2 . (21) 

Their total momentum is given by: 

P = Pi- (22) 

We can observe the same process in the cm. frame, where the two colliding particles have zero total momentum. The 
square s of the 4-momentum in the cm. frame is given by: 

s = PgPo a . (23) 

Since s is Lorentz invariant, it must be: 

s = P a P a = (Ei +m 2 c 2 ) 2 -p\c 2 = (Ei +m 2 c 2 ) 2 - (E\-m\c A ) = m\c A + m 2 2 c A + 2Eim 2 c 2 . (24) 
If the two particles are protons, m\ ~m 2 ~ m p . We then have 

s = 2m p c 2 (Ei + m p c 2 ) . (25) 
The energy of each proton in the cm. frame can be written as: 

E = ym p c 2 , (26) 
in fact, in the cm. frame both of them moves with the same Lorentz factor 7. Since 

s = (2E Q ) 2 , (27) 

we have: 

= y/2m p c 2 (Ei+m p c 2 ) = j Ei+m p c 2 ^ 
2m p c 2 2m p c 2 y 2m p c 2 

An alternative derivation of Eq.(28) starts from the fact that 7 is also the Lorentz factor of the motion of the cm. in 
the laboratory frame, since, we recall, the target proton is at rest in the laboratory frame. The speed v of the cm. in the 
laboratory frame is given by (see e.g. Eq.(l 1.4) in Landau and Lifshitz [189]): 



where we used Eqs.(21)-(22) together with the fact that for two protons m\ = m 2 = m p . From Eq.(29), by definition of 
y, we have: 



(Ei+nipC 2 ) 2 _ J (Ei+nipC 2 ) 2 _ jEi+m p c 2 



l-(v/c) 2 \(E l +m p c 2 f-p\c 2 V E 2 + m 2 c 4 + 2Eim p c 2 -p 2 c 2 ]] 2m p c 2 



where we used the fact that E 2 = p 2 c 2 + m 2 c 4 . As we expected, Eq.(30) is identical to Eq.(28). 

Substituting Eq. (28) in Eq. (20), we get the expression of Fermi for the Lorentz contraction of the volume Q.q in 
the cm. frame with respect to £1: 

/ 2m„c 2 

a = J » n, (3i) 

y Ei + m p c z 

Note that if the energy E\ increase, the volume £2o decreases, as is predicted in special relativity. The parameter volume 
will be therefore energy dependent. Fermi calculated also the cross-section as the area available for collisions around 
the pion cloud 

a to , = nR 2 . (32) 

Substituting X c % of Eq. (19) in Eq. (32), we get a, ot = 6x l(T 26 cm 2 , where h = 1.054 x 10~ 27 erg.s and c = 
2.9979 x 10 10 cm/s a value close to the modern experimental value. 

In treating the collisions of extremely high energy nucleons Fermi make use of thermodynamic laws, instead of 
considering a detailed statistical treatment. The energy density around the colliding nucleons is so high that multiple 
pions as well as antiprotons will be produced. 

From Planck law, the spectral intensity (dimension I v (v,T) — ► dE /dtdAdfldv) of the black body is given as 

Mv,r) = ^^, (33) 

where k is the Boltzmann constant, v the frequency and T the temperature. 

From Stefan-Boltzmann law of the black-body (radiation flux R{T) = aT 4 , dimension dE/dtdA), the energy 
density p(T) (dimension dE/dV) is 



I(T) = ri v (v,T)dv - ^-p(T) = ]-R(T) 
Jo 4n % 

An l"* 4 

p(r) = Ti, I ^ v ^ T ) dv = - c aT > 



(34) 



where n 4 /15 = 6 E 1/n 4 = 6.494 (from Gamma and Riemann Zeta functions, respectively, F(z) and C(s)), O is the 

Stefan-Boltzmann constant and I(T) the spectral intensity integer in all frequency of the black-body. According to 
Fermi: "Consequently the Stefan's law for the pions will be quite similar to the ordinary Stefan's law of the black- 
body radiation. The difference is only in a statistical weight factor. For the photons the statistical wight is the factor, 2, 
because of the two polarization directions. If we assume that the pions have spin zero and differ only by their charge 
±e or 0, their statistical weight will be 3. Consequently, the energy density of the pions will be obtained by multiplying 
the energy density of the ordinary Stefan's law by the factor 3/2." Then, multiplying the energy density (35) by 3/2, 
the energy density via pions is 

PAT) = \ P{T) = l^ {kT)A - (36) 
The total energy of the system (Eq. 25) is divided among pions, protons and anti-protons. Then it is necessary to get 
the energy density via protons and anti-protons. In this case is used Planck law modified (for fermions) as 

2 E 3 

I fermion (E ,T) — J^-J y , (37) 



and ^ 

Pfermion(T) = / If ermion (E,T) — . (38) 

C JO rt 

It is necessary to use the Planck law modified because the protons and anti-protons are fermions, where they obey 
Fermi-Dirac statistical. According to Fermi: "The contribution of the nucleons and anti-nucleons to the energy density 
is given by a similar formula. The differences are that the statistical weight of the nucleons is eight since we have four 
different types of nucleons and anti-nucleons and for each, two spin orientations. A further difference is due to the fact 
that these particles obey the Pauli principle." Then it is necessary to multiply the Stefan-Boltzmann law by fermions 
to 8/2. The process is pp — > n + X, where "X" represent the protons and anti-protons, then the energy density via 
protons and anti-protons is 

(t\ 8 4x5.682 ,^4 
Px{T) = -p fenn ,on{T) = 3 {kT) (39) 

where 6 £ (-l) n+1 /« 4 = 5.682. 

n— 1 

The total energy density p tot of the system during the collision is given by sum p tot = p K + p x , but also it is the 
energy of cm. divided per volume, p tot = ^/s/Clo- Then 

ptot = ^r=Pn + px- (40) 

"0 

Substituting Eqs. (36) and (39) in Eq. (40) 

o ^ _ 3x6.494 4 4x5.682 4 

<* r » 4 = 0152 |?s- < 41 > 

Note that the energy density is frame invariant, p tot = Eq, ,/C1q = E tot /CI, because cancel the Lorentz factors. 

Analogous to the Eqs. (33) and (34) it has the definition of numerical density of the black-body radiation is [see 
190] 

An rUv,T) , %n f°° v 2 dv 
dv = — r ' 



c Jo 



hv c 3 Jo e hv l kT - 1 

n (T) = - f — ) 3 f — = F(3)C(3) ( k Z V 
1 ' n 2 \ch J Jo e*-l n 2 \ch J 



n{T) = 0.243576^-^, (42) 

he 3 



where x = hv /kT. In the case of pions, ( |£§), it is necessary multiply the last equation by 3/2 (analogous Eq. 36), as 

n n {T) = \n(T) =0.365-^, (43) 
2 he 3 

where Fermi [21] got also the last expression. 
Substituting Eq. (41) in Eq. (42), we get 

/ \ 3/4 

nf{^) = 0.0888 [j-L^ . (44) 
We can define the multiplicity of pions per collision (TV™ — > #pions), for high energy, being (using Eq. 31) 

nf = ^^K E = 2 -^Clnf. (45) 



Substituting Ep. (44) in Eq. (45) 

„, _ fm„c 2 Q. LF s\ 1//4 

N H n E (^s) = 0.1777 f p -^- 3 j . (46) 

According Eqs. (18) and (19), 

^ 47Ca 3 /j 3 c 3 

n 'M^7- <47) 



Substituting Eq. (47) in Eq. (46) 



K E (V^) =0.25422 



2 3 

m p c a s 



(nixc 2 ) 3 



1/4 



(48) 



The pionrest masses are different, m n oc 2 = 0.135GeV and m K ±c 2 = 0.139GeV. Fermi got that when considerer the 
conservation of angular momentum, it has the effect of reduction the numbers of pions and nucleons, then he found a 
factor obtained numerically of 0.51. The total energy via pions is divided approximately equal among n°, Tl~ and n + , 
then multiplying and dividing Eq. (48) per, respectively, 0.51 and 3, it gets the n° and ^ ± multiplicities 



N^{Vs) = 0.0432 



2 3 

m n c a s 



1/4 



_(m n0 c 2 y_ 

Doing m^oc 2 = 0.14385m ; ,c 2 and m K ±c 2 = 0.14875m p c 2 in the lest two equations, 



(49) 
(50) 



N»J(VS)= 0.m a y\ U*, (51) 



m„c 



N^(V~s)=0.m a y\h^, (52) 



NZ tal (Vs)=0.546 a y\ (53) 



m n c 



where N^f otal = 2N^± +N™, can note that the lest value is the same in Fermi [21]. 

The center of mass energy is according Eqs.(25)-(27). Substituting in Eq. (48), we get the equation of Fermi to the 
7t° multiplicity for extreme high energies 

N»J(E p ) = OW 4 (l + -^) V4 • (54) 
Doing the same procedure from Eq. (45) to (51), 

N^(VS) = 0.6533 a y ^ /m /f-f 2 . (55) 
71 y/s/nipC 1 

Substituting Eq. (25) in Eq. (55), we get the equation of Fermi to the multiplicity of pions in intermediate energy 
range, 

N$(E P ) = 0.777a 3 / 4 A > . (56) 

y/l +E p /mpC 2 



Modern approach 



Currently the modeling of pp interactions is done through computational codes, Monte Carlo codes such as SIB YLL, 
PHYTHIA, Dpmjet. Kelner [191] presented new parameterizations of energy spectra of secondary particles, % and 77 
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Figure 46. The graphic compares the experimental data (black points) with the SIB YLL code (open points) and with the Kelner's 
expression [191] for the cross section (solid curve). Reproduced from Kelner et al. [191]. 



mesons, gamma rays, electrons, and neutrinos produced in inelastic proton-proton collisions based on the SIBYLL 
code by Lipari [192]. These parameterizations have very good accuracy in the energy range above 100 GeV (see figure 
46 and 47). 

The fit to the pp cross section obtained by Kelner et al. [191] is 



o pp (E p ) = (34.3 + 1.88L + 0.25L 2 ) 



E,h 



mb, 



(57) 



where E p is the incident proton energy in laboratory frame (the same that E p in the last section), L = \n[E p (TeV)} 
and E t h is the minimum threshold energy of the incident proton for production of a pion (E,h = 1.22 GeV) and 
1 mb = 1 mbarn = 10~ 27 cm 2 . 
The function that Kelner got for the multiplicity of pions "N n " is given as 



dN n = F Jt (x,E p )dx, 
F n {x,E p ) = —<j)(x,E p ), 



dx 



1-xPy 



\+kyxPy{\-xb) 



(58) 
(59) 

(60) 



where E K is the total energy via neutral pions productions (secondary pions mesons), x — E n jE p , and B K , j3 r and 
are functions dependent only of E p in the Kelner parametrization [191, see graphic 47]. Figure 47 shows the spectrum 
distribution xF 7t (x,E p ) of 7T° production as a function of x = E n /E p , which is the percentage of incident proton 
energy transferred to the pions. Note in Fig. 47 shoes the small probabilities of the n° productions at small energy 
transformation x < 0.0015 and large energy transformation 1 > x > 0.7 E p . The maximum probability is at x = 0.03, 
indicating the approximate 3% of incident proton energy is transferred via neutral pion n°. Since the experimental data 
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Figure 47. Energy spectra of K and r) mesons from the numerical simulations of the SIBYLL code (histograms) and from the 
analytical presentations given by Eqs. (59) and (60) for energy 0.1 TeV. Reproduced from Kelner et al. [191]. 
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Figure 48. The multiplicity of neutral pions 7t in function of incident proton energy E p for intermediate energy range 50 < E p < 
300 GeV. We compared three results: Kelner-SIBYLL, Fermi approaches in median energy and high energy for a — 5. 



show that the production of neutral pions n is practically the same as the productions of positive n + and negative n~ 
charged pions, then we get that approximate 10% of incident proton energy E p is transferred via pions (E z ). 

Integrating Eq. (58) over x, Kelner obtained [191] the multiplicity (number) of 7T° per pp collision as a function of 

(E p ) =3.92 + 0.83L + 0.075L 2 . (61) 



Comparison between Fermi's and Kelner- SIB YLL's approaches 

We compared the multiplicities of neutral pions production (ft ) obtained the Fermi theoretical approach (Eqs. 54 
and 56) and the Kelner analytical parametrization (Eq. 61) to the SIBYLL code [192]. The analytical parametrization 
of Kelner has good agreement with the SIBYLL code for energy E p > lOOGeV. 
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Figure 49. The multiplicity of neutral pions K° in function of incident proton energy E p for high energy range 300 < E p < 1000 
GeV. We compared two results obtained respectively by Kelner-SIBYLL and Fermi Eq. (54) for a = 5. 



In the figure 48 we compare the Fermi results in high-energy region (Eq. 54) and median-energy region (Eq. 56) 
with the description of Kelner-SIBYLL in intermediate energy range 50 — 300 GeV. We used parameter a = 5, which 
give a very good agreement between Fermi result in high energy and Kelner-SIBYLL result in the range 100 — 300 
GeV. But Fermi result in median energy Eq. (56) gives lower multiplicity of 7T° in this energy range. 

In the figure 49, we use a — 5 and compare the Fermi result for high energy (54) with Kelner-SIBYLL result in the 
energy range 300 — 1000 GeV. We find that Fermi result is about 8% larger than Kelner-SIBYLL result, this indicates 
that the validity of the Fermi approach in extremal high-energy range is in question. Because the cross-section grows 
with E p (see Figure 1), and the a parameter by its definition should become larger with larger as the cross-section 
grows in this energy range. This makes that Fermi result further deviate from Kelner-SIBYLL result. One of reasons 
could probably be the fact that more particles, e.g.,gluons and quarks, are exited and participate the thermalization in 
the Fermi volume ii, as a result, the energy transferred to ^-productions is smaller than that estimated by Fermi with 
three particles proton, neutron and pion. 



Maximum and minimum energy of the Pions 

In this section we will show the energy limits of the pions created via pp interactions and we will apply the limits in 
the work of Blattnig et al. [193]. 

The energy of a pion in the laboratory frame (LF) in function of its energy in center of mass (cm.) frame, is as 

E n = y(Eox + vpoxcos9) , (62) 

where y is the Lorentz factor, "v" the velocity of the pion, Eo K the pion energy in the cm. frame and po n the pion 
momentum in the cm. frame (the index "0" inform cm. frame and without the index "0" give in the LF). We will 
consider c = 1. If cos6 = 1 the pion energy is maximum (E™ ax ), and if cosQ = — 1 the pion energy is minimum (E% ln ). 
Then the maximum and minimum energy of the pion is, 

max 

E™" =y(E 0n ±v PQK ). (63) 
Developing the calculus we obtain the maximum and minimum energy of the pions, 

ET = (l™ P Ep - 2m 2 p + ml +R p ), (64) 

E min = (4^( 2 ">P E p- 2 ">h m l- R p)> if E P> E p (65) 

\m K , if e p <e;, 




Figure 50. a)The confront £ , ™ '"&£ , ™ flA vs E p . b) Application of the minimum (65) and maximum (64) thresholds of the 7l a created 
in the Blattnig parametrization (68). 



where 

R p= W ( E p- m p)l( 2m P E p-2 m p- m x) -16m p m%] 

2m 2 . + 2m„mjr — mi 

K = — ~W7 \ — -~ 1.242 GeV. (67) 

' 2{m p -m K ) 

where E* is the limit (threshold) of energy in the case that E K = m K . 

In Fig. 50a, we can see the behavior of the Eqs. (64) and (65), where they are the maximum and minimum limits of 
the neutral pion energy created. We can note that the minimum limit is approximated of E™" 1 = 0.478 GeV for energy 
of incident proton E„ = 10 4 GeV, where E"™' increase very softly. It is possible note also that the maximum limit tends 
toEJ<=E p -3m p /2. 

Blattnig et al. [193] obtains the spectral distribution of n° [da / dE p (mb /GeV)] in function of the kinetic energy 
of n° created (T n o), for energy of incident proton E p < 50 GeV. Blattnig obtain the spectral distribution and total 
cross section of n° for seven different kinetic energies of incident protons, T p = 0.5, 1.0, 1.9,5.0,9.5,20,50 GeV. The 
analytical function that Blattnig got is 

da ( K 2 K 3 K 4 \ 

dE- p -^{ K ^Tf^W + W)' (68) 

where = —5.8, K% = —1.82, ^3 = 13.5, and K4 = —4.5. The last expression is the analytical fit of a numerical 
integration [193]. In Fig. 50b we apply the maximum (64) and minimum (65) limit energies of the n° in the spectrum 
distribution of Blattnig et al. [193] (Fig. 50b), and we obtained results that contradict their limits of energy, where it 
shows problems in minimum and maximum energies. 

Similar conclusions are being explored using the PYTHIA code [194], and results will be presented soon [195]. 



CONCLUSIONS 

Our current understanding of GRBs is based on a general picture which was presented in a set of letters [4, 5, 6]. On 
that basis a canonical GRB scenario has emerged, with three distinct phases: 

1. The vacuum polarization process occurring in the gravitational collapse to a black hole, and the consequent 
creation of an electron-positron plasma. 

2. An optically thick fireshell characterized by the self-acceleration of such an optically thick electron-positron 
plasma, with the engulfed baryon loading. This phase ends with the reaching of transparency, when the P-GRB 
is emitted. 



3. An optically thin fireshell characterized by an accelerated beam of protons and electrons with a Lorentz y factor 
roughly inversely proportional to the baryon loading, interacting with the CBM. 

This basic scenario is currently evolving in a large number of theoretical details, ranging from a) the thermalization 
process of electrons and positrons after their production, to b) the dynamics of the electron-positron pairs in the 
optically thick phase and their instabilities, to c) the probing and determination of the CBM distribution around the 
gravitationally collapsed object by the interaction of the ultrarelativistic baryons and electrons colliding with the CBM. 

Many of the properties of the observed extended afterglow X- and y-ray emission below ~ 1 MeV has been obtained 
by postulating a thermal spectrum of the emission process in the co-moving frame. This treatment is particularly 
appealing since it allowed to derive explicit analytic formulas to compute all the relativistic transformations between 
the co-moving, the laboratory frame and the arrival time of the signal. This approach has been very satisfactory in 
explaining the overall bolometric luminosity of the sources, its evolution as a function of the arrival time, and the 
CBM filamentary structure. With the improvement of the observational techniques, more and more time resolved 
spectra have been observed and our approach shows some discrepancy in the low and high energy tails of the time 
resolved spectra, although the time integrated spectral distribution is very well recovered. We are currently evolving 
this basic mechanism by assuming some departures from a pure thermal spectral shape in the co-moving frame. 

This analysis has led to a clear identification of sources occurring in CBM with average density n c b m <~ 1 particle/cm 3 
and with n c \, m ~ 10~ 3 particles/cm 3 . The first ones correspond to sources occurring in star forming regions in the host 
galaxy, and the second ones to sources occurring in the galactic halo. This problematic has led to a new understanding 
of the traditional separation between long and short GRBs, as exemplified in these lectures. 

We are now approaching, in view of the new data from the Fermi and AGILE satellites, an analysis of the GRB 
radiation over 1 MeV. It is by now clear that the emission process previously considered of a purely thermal spectrum 
is not appropriate to the description of this high-energy component. 

In parallel, we are currently examining how Fermi ideas [21] have been further developed in large data analysis 
procedures at CERN and other accelerators all over the world (see the last section of this paper). We are going to 
probe, by fitting the observational data, an higher density nodule component in the CBM, which will add to the 
understanding of the CBM itself and will lead very likely to the explanation of GRB sources. 

The major effort now is also directed to the understanding of the process of the electrodynamics of gravitational 
collapse. 
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